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Chapter 2 


Atomic and Nuclear Physics 


1. How many neutrons and protons are there in the nuclei of the following atoms: 
(а) ‘Li, 
[Sol] The nucleus is composed of protons and neutrons and nuclide is expressed by . The total 
number of nucleon, that is protons and neutrons, is equal to A(mass number) - Z(atomic number) 
+ М (neutron number) and atomic number is represented as the proton number. So proton 
number of "Li is 3(atomic number) and the number of neutrons is 4 
(b) "Mg, 
[Sol] 24(mass number) = 12(atomic number) + neutron number. So, the neutron number is 12. 
(c) Xe, 
[Sol] 135(mass number) = 54(atomic number) + neutron number. So, the neutron number is 81. 
(d) Ri, 
[Sol] 209(mass number) = 83(atomic number) + neutron number. So, the neutron number is 126. 


(e) 2? Rn, 


[Sol] 222(mass number) = 86(atomic number) + neutron number. So, the neutron number is 139. 


2. The atomic weight of Со is 58.93319. How many times heavier is "Co than "C? 


[Sol] МС2С) = 12.00000 amu , M(?Co) = 58.93319 amu 
M(?Co) / M(?C) = 58.93319/12.00000 = 4.911 times. 


3. How many atoms аге there in 10g of ^C? 


[Sol] 12g : 0.6022 x 1077 atoms = 10g : X 
X = 102 х 0.6022 x 10™ atoms / 12g = 5.0183 x 10” atoms 


4. Using the data given next and in example 2.2, compute the molecular weights of 
(а) Hogas (6) ЊО (с) H207 


Abundance, а/о 


99.985 1.007825 
0.015 2.01410 





[Sol] (а) Н, gas 
M(H, gas) = { (99.985 x 1.007825)/100 + (0.015 х 2.01410)/100} x 2 
= 2.015952 amu 


(b) Н.О 

М(Н,О) = М(Н, + M(O) 
2.015952 +0.01(99.759 х 15.99492 + 0.037 x 16.99913 + 0.204 x 17.99916) 
= 18.05532 amu 


(c) H20, 
М(Н,О,) -М(Н,) + 2M(O) 
= 2.015952 + 31.99876 
= 34.014712 amu 


5. When Н» gas is formed from naturally occurring hydrogen, what percentages of the molecules have molecular 
weights of approximately 2, 3, and 4? 


[Sol] (a) In the case of being molecular weights of 2 
'H +'Н > Н, 


0.99985 x 0.99985 = 0.9991 = 99.91 % 


(b) In the case of being molecular weights of 3 
'H +7H > Н, 


0.99985 x 0.00015 = 0.000149977 = 0.0149977% 
(c) In the case of being molecular weights of4 


0.99985 x 0 + 0.00015 x 0.00015 = 0.00000002 % 


6. Natural uranium is composed of three isotopes, ^^U, 7U, ?*U, Their abundances and atomic weights are 
given in the following table. Compute the atomic weight of natural uranium. 


Abundance, 96 
0.0057 234.0409 


235.0439 
99.27 238.0508 





[Sol] M(U) = 0.01(0.0057 x 234.0409 + 0.72 x 235.0439 + 99.27 x 238.0508) 
= 238.0187 amu. 


7. A beaker contains 50g of ordinary (1.е., naturally occurring) water. 
(a) How many moles of water are present? 

(b) How many hydrogen atoms? 

(c) How many deuterium atoms? 


[Sol] (a) М(Н,О) = 18.015g/mole 
509/18.0015g-mole' = 2.775mole 


(b) 2.775mole x 0.6022 x 102 atoms/mole х 2 
= 3.3422 x 10” 


(c) Abundance of Н is about 0.0151 % 
3.3422 x 107" x 0.01 x 0.0151 25.046 х 1077 atoms 


9. The glass in Example 2.1 has and inside diameter of 7.5 cm. How high does the water stand in the glass? 


[Sol] N(H) = 6.6 x 10" atoms 
N(H;O) = 0.5 N(H) 
= 3.3 x 10" atoms 
3.3 x 10" atoms / (0.6022 x 10" atoms: mole” ) = 5.4799 mole 
p(density) = m(mass)/v(volume) 
m= 5.4799 mole x 18g/mole = 98.6382g 
р(Н,О) 298.6382g/ (3.75? x л x Н ст?) 
Н = 2.2327 cm 


9. Compute the mass of a proton in amu. 


[Sol] The mass of proton is 1.67261 x 10-^*g 


Mp = 1.67261 x 10g x (amu/1.66053 х 10-”g) 
= 1.007277 amu 


10. Calculate the mass of a neutral atom of 21) 
(a) in amu 
(b) in grams 


[Sol] M(^*U)2235.0439 amu = 235.0439 amu x 1.66053 x 107*g/amu 
= 3.903 х 10726 


11. Show that 1 amu is numerically equal to the reciprocal of NA 


[Sol] 1 ати = 1/12 x m(12C) 
=1/12 x 1.99264 х 1026 
= 1.66053 х 10 g 
= 1/(1/1.66053 х 1077 ) р 
= 1 / 6.02216 x10? е 
= 1/ №о 


12. Using Eq. (2.3), estimate the radius of the nucleus of ^*U. Roughly what fraction of the ^*U atom is taken 
by the nucleus? 


[sol] Eq. (2.3) В, =1.25fm x А!” = 1.25 x 10? x 2381” = 7.7464 x 1077 ст 
R;(Radius of atom) = 2 x 10? cm 


13. Using Eq. (2.3), estimate the density of nuclear matter in g/cm”; in Kg/m'. Take the mass of each nucleon to 
be approximately 1.5 10% g. 


[sol] Nuclear matter = ^X 
R(Radius of nucleus) = 1.25fm x А' = 1.25 x 19? x A'? 


14. The planet earth has a mass of approximately 6 10“ kg. If the density of the earth were equal to that of 
nuclei, how big would the earth be? 


[sol] Density = 1833.5x10'* kg/m? (2.13 ans) 


R = 198.42m (Radius of earth = 6400km) 


15. The compete combustion of 1 kg of bituminous coal releases about 3 10’ J in heat energy. The 
conversion of 1 g of mass into energy is equivalent to the burning of how much coal? 


[sol] E. = Energy of 1 kg of bituminous coal = 3x10’ J 


16. The fission of the nucleus of ^?U releases approximately 200 MeV. How much energy (in kilowatt-hours 


and megawatt-days) is released when 1 g of ^?U undergoes fission? 


[sol] № = 0.6022045x10?^/mol 


17. Compute the neutron-proton mass difference in MeV. 
[sol] Neutron : 
Proton : 


Neutron-proton mass difference = 1.28MeV 


18. An electron starting from rest is accelerated across a potential difference of 5 million volts. 
(a) What 18 its final kinetic energy? 


[sol] 
(b) What is its total energy? 


[sol] 


(c) What is its final mass? 


[sol] 


19. Derive Eq. (2.18). 


[sol] Eq. (2.5) 





Eq.(2.18) 


20. Show that the speed of any particle, relativistic or nonrelativistic, is given by the following formula: 


— — , where Erest and Etotal are its rest-mass energy and total energy, respectively, and c is the 


speed of light. 


[sol] 


2]. Using the result derived in Problem 2.20, calculate the speed of a 1-MeV electron, one with a kinetic energy 
of 1 MeV. 


[sol] 


22. Compute the wavelengths of 1-MeV 
(a) photon 


[sol] 
(b) neutron. 


[sol] 


23. Show that the wavelength of a relativistic particle is given by 


? 


where is called the Compton wavelength. 
[sol] When body is in motion, it's mass increase relative to an observer : ————----- (Eq 2.5), 
So 


Relativistic momentum : p = mv = ————-—--(Eq.2.15) 


Som ———, —— 


Rearrange this relation. 


so pc-2  --- 
р = -, replace p in eq -- (1) => е = 


60 = ———-_ — 


24. Using ће formula obtained in Problem 2.23, compute ће wavelength of a 1-MeV electron. 


[sol] in prob 2.23, — > 


[kg] (2.9979 [] $81.87 [joule] 
lev = 1.60219 joule, 
81.87 1.60219 510.9950 0.511[MeV] 


ен 


25. An electron moves with a kinetic energy equal to its rest-mass energy. Calculate the electron's 
(a) total energy in units of 


[sol] = - , so =2 


(b) mass in units of 


[501] m=, m=2 


(c) speed in units of c 


[sol] = ----- Eq 2.8 








=> = 1/4 





m=2,so 2= 


—= —2.9979 


(d) wavelength in units of the Compton wavelength 


26. According to Eq. (2.20), a photon carries momentum, thus a free atom or nucleus recoils when it emits a 
photon. The energy of the photon 18 therefore actually less than the available transition energy(energy between 
states) by an amount equal to the recoil energy of the radiating system. 

(a) Given that E is the energy between two states апа .1$ the energy of the emitted photon, show that 


Where M 18 the mass of the atom or nucleus 


[sol] The momentum of zero mass particle is 


р = E/c ----(Е 4.2.20) 
The momentum of the particle having nonrelativistic energy is 
р= ----(Eq.2.15) 


(b) Compute E - for the transitions from the first excited state of atomic hydrogen at 10.19eV to ground 
and the first excited state of at 4.43 MeV to ground (see Figs. 2.2 and 2.3) 


[sol] i) For Hydrogen 
Е- = —— 


ii) For Carbon 


27. The first three excited states of the nucleus of are at 0.158 MeV, 0.208MeV, and 0.403 MeV above 


the ground state. If all transitions between these states and ground occurred, what energy -rays would be 
observed? 


[sol] 
0.403MeV Possible observed -rays are .. 





0.208MeV i) 
0.158MeV ii) 
Ground iii) 

iv) 

v) 

vi) 


28. Using the chart of the nuclides, complete the following reactions. If a daughter nucleus is also radioactive, 
indicate the complete decay chain. 


(a) 

[sol] (stable) 

(b) 

[sol] 

(c) 

[sol] 

(d) 

[sol] 
29. Trittum( decays by negative beta decay with a half-life of 12.26 years. The atomic weight of is 
3.016 

(a) To what nucleus does decay? 

[sol] decay: neutron proton + -anti-neutrino so it will be cause increase in proton number 


and same atomic mass So decay of ‚ the proton number increase +1 so it will be 2, and same 
atomic mass 3, 


(b) What is the mass in grams of 1mCi of tritium? 


[sol] ----(Ед.2.23), using Eq.2.26 —— 


111(1 = 3.7 disintegrations / second = 
Total number of ,n(t) 
n(t) 23.7 2.06 


mass = —= 


30. Approximately what mass of ( has Фе same activity as 1g of 


[sol] Activity of 1g of 
= —p———— 
= 4.19 | 
some mass of activity same as 1g of Co activity 
Е 
9 = -5,50--- 
— бо,Х= --8:2197990088 


31.Carbon tetrachloride labeled with 15 sold commercially with an activity of 10 millicuries рег 
millimole(10mCi/mM). What fraction of the carbon atoms is ? 
[sol] Activity = 10mCi/mM = 3.7 =n(t) :— is 5730yr 

n(t)2 ——— — — ——7 9.6458 


32. Tritiated water (ordinary water containing some for biological applications can be purchased in 
1- ampoules having an activity of 5mCi per . How many atoms per . What fraction of the water 
molecules contains an atom? 


[sol] , atomic mass : 20.04g/mole 
volume of 1mole of 
———— — —= 3.6445 disintegrations per second 


9 


n(t) = 


33. After the initial cleanup effort at Three Mile Island, approximately 400,000 gallons of radioactive water 
remained in the basement of the containment building of the Three Mile Island Unit 2 nuclear plant. The 
principal sources of this radioactivity were at 156 Ci/ and at 26 Ci/. How many atoms 

per of these radionuclides were in the water at that time? 


[sol] For 


n(t)= 


For 


n(t)- 
number of total atoms at that time per 
of atoms and of atoms 


34. One gram of Ra is placed in a sealed, evacuated capsule 1.2cm? in volume. 
(a) At what rate does the helium pressure increase in the capsule, assuming all of the a-particles are 
neutralized and retained in the free volume of the capsule? 


[sol] pressure rate — activity 
(per volume), — —— 


(b) What is the pressure 10 years after the capsule is sealed? 


[sol] 


35. Polonium-210 decays to the ground state of *°°Pb by the emission of a 5.305-MeV a-particle with a half-life 
of 138 days. What mass of 210Po is required to produce IMW of thermal energy from its radioactive decay? 


[sol] 


=0.0314056W/Ci 
To produce 1MW(thermal), 


| — We need this reactivity . 


36. The radioisotope generator SNAP-9 was fueled with 475 g of ^*PuC (plutonium-238 carbide), which has а 
density of 12.5 g/cm’. The Ри has a half-life of 89 years and emits 5.6MeV per disintegration, all of which 
may be assumed to be absorbed in the generator. The thermal to electrical efficiency of the system is 5.4%. 
Calculate 
(a) the fuel efficiency in curies per watt (thermal); 


[sol] 
(b) the specific power in watts (thermal) per gram of fuel; 


[sol] —— — — — — M — 
= 4.763 


(с) the power density in watts (thermal) per ст”; 


[sol] —— — — ——— watts/ 


(d) the total electrical power of the generator. 


[sol] 


37. Since the half-life of 270 (7.13 x10* years) is less than that of 2380 (4.51x10? years), the isotopic 
abundance of 270 has been steadily decreasing since the earth was formed about 4.5 billion years ago. How 
long ago was the isotopic abundance of ^?U equal to 3.0 а/о, the enrichment of the uranium used in many 
nuclear power plants? 


[sol] 


(Primary) Abundance of 





9 


——0.0072— — — 


38. The radioactive isotope Y is produced at the rate of R atoms/sec by neutron bombardment of X according to 
the reaction 
X(n, y) Y. 
If the neutron bombardment is carried out for a time equal to the half-life of Y, what fraction of the 
saturation activity of Y will be obtained assuming that there is no Y present at the start of the bombardment? 


[sol] = 
Assuming that there is no Y present at the start of the bombardment, 


If the neutron bombardment is carried out for a time equal to the half-life of Y, 


= 0.5В 


39. Consider the chain decay 
A—> B> С >, 
with no atoms of В present at t=0. 
(a) Show that the activity of B rises to a maximum value at the time tm given by 


? 


at which time the activities of A and B are equal. 


[sol] ——— 


— -ð 


(b) Show that, for t<tm, the activity of В is less than that of A, whereas the reverse is the case for t> ty. 


[sol] 
— Ф : 
< 090: — 
If 
If 
If 
If 


40. Show that if the half-life of B is much shorter than the half-life of A, then the activities of A and B in 
Problem 2.39 eventually approach the same value. In this case, A and B are said to be in secular equilibrium. 


[$01] _ — — 





If 


41. Show that the abundance of ^^U can be explained by assuming that this isotope originates solely from the 
decay of ^*U. 


[sol] 


SO -——> 


24.1days 6.75years 
— Halflife of 15 much longer than others (Th, Pa). 


So depends on decay of 


42. Radon-222, a highly radioactive gas with a half-life of 3.8 days that originates in the decay of ^"U (see chart 
of nuclides), may be present in uranium mines in dangerous concentrations if the mines are not properly 
ventilated. Calculate the activity of Ви іп Bq per metric ton of natural uranium. 


[sol] ——— 














43. According to U.S. Nuclear Regulatory Commission regulations, the maximum permissible concentration of 
radon-222 in air in equilibrium with its short-lived daughters is 3 PCi/liter for nonoccupational exposure. 


This corresponds to how many atoms of radon-222 per сш"? 


[sol] — = 





Short-lived daughters 


= 3min + 27min + 20min + 180microsec 
= 3000.18s 


44. Consider again the decay chain in Problem 2.39 in which the nuclide A is produced at the constant rate of R 
atoms/sec. Derive an expression for the activity of B as a function of time. 


[sol] 1. — -ANA +R= -(An, = К) 


If t = 0, ny = no 
П = ct - 
C+R=A) по 
С =Лш-В 


= По - — + ын 
=m-  &-) 

2.A5 BC 
Е Agng = МАП А 
ПА = Nao - a-) 
—— рр = Ал Nao - R(1 - ) шинйшэндиливнэийн а 
2-1. homogeneous 
--Т Agng =0 
dng = -)gnpd, 

в = d, 

= -Apt + C 
Пв = C 
if t — 0, Пв = Ngo 
Ngo = С 
Np = Пво 
2-2. non-homogeneous 
Let, Apt =I 
)gdt = dI 
n А 
Пв = Пво 
Пв = Пво 
(ав) = п'в + ng = 

= n'g + NgÀg 

Eq. a x 
— = Agng = АА Nao - R(1 - ) = ng) 
Пв = dt - dt 4C 
Пв = dt - dt +C 


= dt- +С 
-—R — уже 
ng = — = re 
if t= 0, ng = 0 
—RC — }+E=0 
C=R- }— —— 
ng = —R¢ - HR - э— — 
= (+ Өс €-1) 
3. homogeneous + non-homogeneous 
Пв = Пво + + + є-1) 
activity Gg = Ово + {=} + + €-1) 


45. Complete the following reactions and calculate their О values. [Note: The atomic weight of “С is 
14.003242.] 
[a] "He(p, d) 


[sol] “Не(р, d) = *He(p, d) He : “Не + ЇН 5 *He +?Н 
О={[М( Не) + M(H)] - [МС Не) + MÓH)]) х 931.502 MeV «Eq. (2.38)> 
= [(4.002603 + 1.007825) — (3.016029 + 2.014102)] x 931.502 MeV 
= -18.3534 MeV 


[b] “Ве(о, n) 


[sol] “Ве(о, п) = Ве(а, п)! °С: Ве + “He > "C + п 
О = [МСВе) + МСНе)] — [МС?С) + M(n)]} x 931.502 MeV 
= [(9.012182 + 4.002603) — (12 + 1.008665)] х 931.502 MeV 
= 5.70079 MeV 
[c] "N(n, p) 


[sol] "N(n, р) = "N(n, p)“C : ИМ + n o "C +'H 
О ={[MC N) + M(n)] - IM(*C) + МО Н)]} x 931.502 MeV 
= [(14.003074 + 1.008665) — (14.003242 + 1.007825)] x 931.502 MeV 
= 0.625969 MeV 
[d] “шд, p) 


[sol] (а, p) = Га, р)" : а + 2Н — 1а + ‘Н 
О = ([[M(?In) + MÓH)] - [МОТ + M(!H)] x 931.502 MeV 
= [(114.903879 + 2.014102) — (115.905262 + 1.007825)] x 931.502 MeV 
= 4.55877 MeV 
[e] ““РЬСу, n) 


[sol] ""Pb(y, n) = ?"pp(y, ny 6pp : "pp + ү — Pb + п 
О = ([M("Pb) + М(у)] - I[MC' Pb) + M(n)]) x 931.502 MeV 
= [(206.975880 + 0) — (205.974449 + 1.008665)] x 931.502 MeV 
= -6.73849 MeV 


46. (a) Compute the recoil energy of the residual, daughter nucleus following the emission of a 4.782-MeV a- 
particle by бра. 


[sol] "*Ra > ??Rn + a : "Ra — Rn + “Не 
Recoiling Energy - Q(E) — E, 
О = [MC Ra)] - [MC”Rn) + МСНе)]} х 931.502 MeV 
= [226.025402 — (222.017570 + 4.002603)] x 931.502 MeV 
= 4.87082 MeV = 4.871 MeV 
Recoiling Energy = E — E, = (4.871 — 4.782) MeV = 0.089 MeV 
(b) What is the total disintegration energy for this decay process? 


[sol] Total disintegration energy = Q value 
= 0.089 MeV 


47.In some tabulations, atomic masses are given in terms of the mass excess rather than as atomic masses. The 
mass excess, A, is the difference 
A=M-A, 
Where М is the atomic mass and A 15 the atomic mass number. For convenience, A, which may be positive or 
negative, is usually given in units of MeV. Show that the Q 
О = (A, + Ay) – (Ac + Aa), 


[sol] О = М, + My) - (M, + Mg)]931 MeV «Eq. (2.38)> 
= (А, + Ag + Ap + Ap) — (Ас +А. + Aa + Ag)] <А = M — A and usually given in units of MeV» 
= (Aa + Ay) — (Ac + Aa) <А, +Аъ = Ae + Aa? 


48. According to the tables of Lederer et al. (See References), the mass excesses for the (neutral) atoms in the 
reaction in Example 2.8 are as follows: ACH) = 14.95 MeV, ACH) = 13.14 MeV, A(n) = 8.07 MeV, and A(*He) = 
2.42 MeV. Calculate the Q value of this reaction using the results of Problem 2.47. 


[sol] a- b— с+а 
ЗН +7H — Не + n 
О = (Aa + Ab) - (Ac + Ad) «Problem 47» 
= (АЁН) + ACH) - (A Не) + A(n))] 
= [(14.95 + 13.13) — (2.42 + 8.07)] MeV 


= 17.59 MeV 


49. The atomic weight of 20РЬ is 205.9745. Using the data in Problem 2.35, calculate the atomic weight of "Ро. 
[ Caution: See Problem 2.46] 


[sol] E, 2 5.305 MeV «Problem 2.35» 
210р), вэ" 206pp +a: 210Pþ = 206pp + ^H e 
0-1МС Pb)] - [МСУРЬ) + MCHe)]) x 931.502 MeV 
= [MP Pb) - (205.9745 + 4.002603)] x 931.502 MeV = 5.305 MeV <Recoiling Energy = 0» 
M(?*Pp) = 209.9828 


50. Tritium (CH) can be produced through the absorption of low-energy neutrons by deutrerium. The reaction is 
^H en > Н + y, 
Where the y-ray has an energy of 6.256 MeV. 
(a) Show that the recoil energy of the ^H nucleus is approximately 7 KeV. 


[sol] E, 2 6.256 MeV 


Recoiling Energy = МУ? = —<р=шу> = -«p- > 
= (6.256 MeV)? / (2 x МСН) x 931.502 MeV) 
= (6.256 MeV)? / (2 х 3.016049 х 931.502 MeV) 
= 0.006965 MeV = 7 KeV 
(b) What is the Q value of the reaction? 


[sol] Q -1МСН) + M(n)] - МСН)}х 931.502 MeV «Eq. (2.38)> 
= [(2.014102 + 1.008665) — 3.016049] x 931.502 MeV 
= 6.25783 MeV 
(c) Calculate the separation energy of the last neutron in ^H. 


[sol] Separation Energy - Binding Energy of the last neutron 
E, = [M, + M(*!Z) - M(“Z)] x 931.502 MeV «Eq. 2.45» 
= [M, + МСН) - MCH)] x 931.502 MeV 
= 6.25783 MeV 
(d) Using the binding energy for ^H of 2.23 MeV and the result from part (c), compute the total binding 
energy of ^H. 


[sol] ВЕс z; = [ZM(H) + NM, - M(*Z)] х 931.502 MeV 
ВЕ 1.2) = [ZM(H) + (N-1)M, - M(*"Z)] x 931.502 MeV 
BE, 7; - ВЕ дл) =[ZM(H) + NM, - M(^Z) - ZM(H) - (N-1)M, + M(*'Z)] x 931.502 MeV 
= [M, + M(*"Z) - M(à^Z)] x 931.502 MeV = E, 
BE, 7; = E, + BE44,7, = BECH) = Е, + ВЕСН) = (6.25783 + 2.23) MeV = 8.48783 MeV 


51. Consider the reaction 
?Li(a, р) Be 
Using atomic mass data, compute: 
(a) the total binding energy of °Li, "Be, and “He: 


[sol] ВЕ(а) = [Z,M(H) + NM, - M,] x 931.502 MeV 
M(H) = 1.007825 
M, = 1.008665 
ВЕСТА) = (3 x 1.007825 + 3 x 1.008665 — 6.015122) x 931.502 MeV 


= 31.9952 MeV 
BE( Be) = (4 x 1.007825 + 5 x 1.008665 — 9.012182) x 931.502 MeV 


= 58.1658 MeV 
ВЕСНе) = (2 x 1.007825 + 2 x 1.008665 — 4.002603) x 931.502 MeV 


= 28.2962 MeV 
(b) the Q value of the reaction using the results of part (a). 


[sol] Q = [(M, + My) - (М, + Mj)]c? XEq. 2.37» 


МСТ + MC He) - МСВе) - Mplc? 
[3M(H) + 3Mn - BE(Li)/c? + 2МСН) + 2Mn - BE(He)/c? 








-BE(Be)-BE(Li)-BE(He) 
= (58.1658 — 31.9952 — 28.2962) MeV 
= -2.1256 MeV 


52. Using atomic mass data, compute the average binding energy per nucleon of the following nuclei: 
(а) Н 


[sol] ВЕСН) = (1 x 1.007825 + 1 х 1.008665 — 2.014102) х 931.502 MeV 
= 2.22443 MeV 
ВЕСН)/2 = 1.11221 MeV 
(b) “He 


[sol] ВЕСНе) = (2 x 1.007825 + 2 x 1.008665 — 4.002603) x 931.502 MeV 
= 28.2962 MeV 
ВЕСНе)/4 = 7.07406 MeV 
(c) 120 


[sol] ВЕС2С) = (6 x 1.007825 + 6 x 1.008665 — 12) x 931.502 MeV 
= 92.1628 MeV 
BE(’C)/12 = 7.68023 MeV 
(d) "V 


[sol] BEC! V) = (23 х 1.007825 + 28 x 1.008665 — 50.943964) x 931.502 MeV 
= 445.846 MeV 
BE(*!V)/51 = 8.74207 MeV 
(e) 138p4 


[sol] BE( Ва) = (56 x 1.007825 + 82 х 1.008665 — 137.905242) x 931.502 MeV 
= 1158.31 MeV 
BE(5Ba)/138 = 8.39356 MeV 
(f) 2551) 


[sol] BEC°U) = (92 х 1.007825 + 143 x 1.008665 — 235.043922) х 931.502 MeV 
- 1783.89 MeV 
BE(^?^Uy/235 = 7.59104 MeV 


53. Using the mass formula, compute the binding energy per nucleon for the nuclei in Problem 2.52. Compare 
the results with those obtained in that problem. 
(а) Н 


[sol] М = NM, + ZM, - aA + BA”? + yZ?/A'? + (А – 27)2/А <Ед. 2.50» 
M = NM, + ZM, - aA + BA?” + yZ^/A'? + à(A – 27)?/А + б XEq. 2.51» (Nuclei with odd numbers 
of neutrons and odd numbers of protons would thus be less strongly bound together.) 
BE = aA - ВА? - у22/А! - б(А – 22)/A – б 
BECH) = (15.56 x 2 — 17.23 х 22° — 0.697 x 12/213 — 23.285 x (2 — 2 x 1)?/2) MeV <N=1, Z=1> 
= 3.21587 MeV 
ВЕ(?Н)/2 = 1.60794 MeV (> 1.11221 MeV) 
(b) “Не 


[sol] ВЕ( Не) = (15.56 x 4 — 17.23 x 47° — 0.697 x 27/4? — 23.285 x (4 - 2 x 2)?/4 — 12.0 x -1) MeV <N=2, 
Z=2> 
= 29.0668 MeV 
BE(*He)/4 = 7.2667 MeV (> 7.07406 MeV) 


[sol] BE(?C) = (15.56 х 12 — 17.23 x 122° — 0.697 x 67/12'% — 23.285 x (12 — 2 х 6)2/12 — 12.0 x -1) MeV 


«N26, 7-6» 
= 97.4493 MeV 
BE("C)/12 = 8.12078 MeV (> 7.68023 MeV) 
(d) Sly 


[sol] ВЕСУ) = (15.56 х 51 — 17.23 x 512° — 0.697 х 237/51! — 23.285 x (51 — 2 x 23)2/51 — 12.0 x 0) 
MeV <N=23, Z=28> 
= 445.765 MeV 
BE(*!V)/51 = 8.7405 MeV (< 8.74207 MeV) 


(е) “Ва 


[sol] BE( Ва) = (15.56 x 138 — 17.23 х 1382? — 0.697 x 567/138" — 23.285 x (138 — 2 x 56)2/138 — 12.0 х 
-1) MeV <N=56, 7-82» 
- 1162.12 MeV 
ВЕС”Ва)/138 = 8.42117 MeV (> 8.39356 MeV) 


(f) 2551 


[sol] ВЕ(20) = (15.56 x 235 — 17.23 х 235°” — 0.697 х 922/235! — 23.285 x (235 — 2 x 92)2/235 — 12.0 х 
0) MeV <N=92, Z=143> 
= 1786.75 MeV 
BE(^?Uy/235 = 7.60319 MeV (>7.59104 MeV) 


54. Compute the separation energies of the last neutron in the following nuclei: 
E,-[M, + M(*!Z) - M(^Z)]931 Mev/amu 
(a) *He 
[sol] E,=[1.008665 + 3.016029 — 4.002603]931=20.567 MeV 


(b) 14 
[sol] Е.=[1.008665 + 6.015122 — 7.016004]931=7.246 MeV 


(c) UO 
[sol] E,=[1.008665 + 15.99491 — 16.99913]931=4.141 MeV 


(d) v 
[sol] E,=[1.008665 + 49.947163 — 50.943964]931=11.045 MeV 


(e) 208pp 
[sol] E,=[1.008665 + 206.97588 — 207.97664]931=7.3633 MeV 


[sol] E,=[1.008665 + 234.040945 — 235.043922]931=5.296 MeV 


55. Derive Eq.(2.53).[Hint: Try taking the logarithm of Eq.(2.52) before differentiating. | 


[sol] —— 














56. What is 1 atmosphere pressure in units of eV/cm ? [Hint: At standard temperature and pressure (0С and 1 
atm), 1 mole of gas occupies 22.4 liters.] 


[sol] The Eq.(2.58), P=NkT 


57. Calculate the atom density of graphite having density of 1.60 g/cm’. 


[sol] In Eq.(2.59) —-50 


58. Calculate the activity of 1 gram of natural uranium. 


[sol] 
Abundance, а/о Half life (T 1/2) 
U-234 0.0057 234.0409 2.4610E5 y 


U-235 235.0439 7.03810Е8 y 
U-238 99.27 238.0508 4.6810Е9 у 


Atomic weight of natural uranium is; 





Total number of natural uranium nuclei in 1 gram of uranium is; 


In Eq. (2.26), A=In2/T 12 


59. What is the atom density of ^?U in uranium enriched to 2.5 а/о in this isotope if the physical density of the 
uranium is 19.0 g/cm?? 


[sol] 


60.Plutonium-239 undergoes a-decay with a half-life of 24,000 years. Compute the activity of 1 gram of 
plutonium dioxide, **’PuO. Express the activity in terms of Ci and Ва. 


[sol] 





The weight percent of plutonium is 


The weight percent is almost the same with а/о. 


61. it has been proposed to use uranium carbide (UC) for the initial fuel in certain types of breeder reactors, with 
the uranium enriched to 25 w/o. The density of UC is 13.6 g/cm”. 
(a) What is the atomic weight of the uranium? 


[sol] — — — 


(b) What is the atom density of the гэр) 


[sol] Molecular weight of UC is; 


62. Compute the atom densities of ^?U and ^*U in ОО, of physical density 10.8 g/cm?’ if the uranium is 
enriched to 3.5 w/o in ^?U. 


[sol] Eq.(2.65) - — — 


ңе i — 


The weight percent of U in UO; is; 


63. The fuel for a certain breeder reactor consists of pellets composed of mixed oxides, UO, апа РиОЬ:, with the 
PuO, comprising approximately 30 w/o of the mixture. The uranium is essentially all ^*U, whereas the 
plutonium contains the following isotopes: Ра (70.5 w/o), ““Ри(21.3 w/o), ^''Pu(5.5 w/o), and ^"Pu(2.7 w/o). 
Calculate the number of atoms of each isotope per gram of the fuel. 


[sol] 





Chapter 3 


Interaction of Radiation 
with Matter 


1. Two beams of 1-eV neutrons intersect at an angle of 90°. The density of neutrons in both beams 152 10* 
neutrons/cm 


(a) Calculate the intensity of each beam. 


[sol] = 1.383 10° cm/s «Eq. 2.12» 
= 1.383 10° cm/s 
= 1.383 10*cm/s 
О = nv XEq. 3.20» 
=(2 105neutrons/cm^) (1.383 10*cm/s) 
= 2.766 10" neutrons/cm?-sec 


(b) What is the neutron flux where the two beams intersect? 


[sol] 


the neutron flux in the region where they intersect «law of cosines> 


= 3.289 10" neutrons/cm?-sec 


2. Two monoenergetic neutron beams of intensities I; = 2 1079 neutrons/cm^-sec and b = 1 10° 


neutrons/cm"-sec intersect at an angle of 30°. Calculate the neutron flux and current in the region where they 
intersect. 


[sol] 


ca \ 


the neutron flux in the region where they intersect «law of cosines> 


= 0.969 10” neutrons/cm?-sec 


: 1 2 : : 2. 
3. A monoenergetic beam of neutrons, Ø =4 10? neutrons/cm?-sec, impinges on a target 1 cm in area and 


0.1 cm thick. There are 0.048 10% atoms per ст” in the target, and the total cross-section at the energy of the 
beam is 4.5 b 


(a) What is the macroscopic total cross-section? 
[sol] Уч =N Ot 

= (0.048 10?) (4.5 102) 

= 0.216 cm” 


(b) How many neutron interactions per second occur in the target? 


[sol] interactions rate = o,INAX (1: intensity of the beam, A : area, X : thickness) «Eq. 3.2» 


= НАХ 
-0.216 ст! (4 10 neutrons/cm?-sec) 1cm? 0.1сш 


=8.64 10° interactions(interaction neutrons)/sec 


(c) What 18 the collision density? 


[sol] collision density Е = INo; (Eq. 3.7» 
= ІУ, 
= (4 10" neutrons/cm?-sec) 0.216 cm” 
= 8.64 10° collisions(collision neutrons)/cm*-sec 


4. The B-emitter ^AI (half-life 2.30 тіп) can be produced by the radiative capture of neutrons by *’Al. The 


0.0253-eV cross-section for this reaction 18 0.23 b. Suppose that a small, 0.01-g aluminum target is placed in a 
beam of 0.0253-eV neutrons, @ = 3 10° neutrons/cm7-sec, which strikes the entire target. Calculate 


(a) the neutron density in the beam; 


[sol] I = nv «Eq. 3.1» 


n=I/v 
=(3 10° neutrons/cm7-sec) / (2.2 107 cm/sec) 


<the speed of thermal neutron = 2.2 km/sec) 
= 1363.64 neutrons/cm? 


(b) the rate at which ^AI is produced; 


[sol] Production rate = INo, 
N = pN,/M «Eq. 2.59» 
INo, = (3 10*neutrons/cm?-sec) {(0.01 g/em3 0.6022 102)/26.9815386) 
0.23 10” ст? 
«M» = 26.981538> 
=1.54 10! atoms/cm*-sec 


(c) the maximum activity (in curies) that can be produced in this experiment. 





[sol] = INC"ADo, — МСА!) 
Where IN("ADo, = С and NAI = N, = C- AN € 
—= С -AN 
—+АМ = С 
€N) = С 


- N(0) = €- € 
- 6-4) 

= @-) 

a(t) =С(1-) 


OMaximum activity = 
= C = IN ADoc; 

= 1.54 10“ atoms/cm*-sec 

= 4.16 107” С! <1Сі = 3.7 10" disintegrations/sec> 


5. Calculate the mean free path of 1-eV neutrons in graphite. The total cross-section of carbon at this energy is 


4.80 


[sol] >, =N с, 
= (0.08023 107) (4.8 107^ «the density of carbon > Table П.3» 
= 0.3851 cm” 
mean free path À = 1 / X, = 2.5967 cm 


б. A beam of 2-MeV neutrons is incident on a slab of heavy water (D20). The total cross-sections of deuterium 
and oxygen at this energy are 2.6 b and 1.6 b, respectively. 
(a) What 18 the macroscopic total cross-section of D;O at 2MeV? 


[sol] X, 2 N 6o, 
N(D,O) = 0.03323 10” molecules/cm? «Table П.З» 
o of D-O = (2 o of deuterium) + с of oxygen 
=(2 2.6)+1.6=6.8b 
У,-(0.03323 107) (6.8 107) = 0.225964 ст" 


(b) How thick must the slab be to reduce the intensity of the uncollided beam by a factor of 10? 


[sol] I = I? «Eq. 3.11» 


- 2.x = -|n10 
x = In10 / X, = 1010 / 0.225964 cm” = 10.19 cm 


(c) If an incident neutron has a collision in the slab, what is the relative probability that it collides with 
deuterium? 


[sol] relative probability = X, of deuterium / X, of heavy water 
={N(D,O) (2 соғ deuterium)} / 0.225964 
= {(0.03323 107) (2 2.6 107%4)}/ 0.225964 
= 0.7647 


7. A beam of neutrons is incident from the left on a target that extends from x = 0 to x =a. Derive an expression 
for the probability that a neutron in the beam will have its first collision in the second half of the target that is in 
the region a/2 < x <a. 


[sol] р(х) 4х = dx <page 59> 
P = = T 


! 
+ 


9. Boral 18 commercial shielding material consisting of approximately equal parts by weight of boron carbide 
(В.С) and aluminum compressed to about 95% theoretical density (2.608 g/cm’) and clad with thin sheets of 
aluminum 0.25 cm thick. The manufacturer specifies that there are 0.333g of boron per cm’ of а boral sheet 
0.457 cm in overall thickness. What is the probability that a 0.0253-eV neutron incident normally on such a 
sheet will succeed in penetrating it? 


[sol] 


Penetrating probability 
- the probability that a neutron can move through this distance without having a collision 
= «page 58» 
M(B4C) 24 10.811 + 12 = 55.244 
mass of C per cm? = —9:333g = 0.07232 
mass of В,С per ст? = 0.333g + 0.0723g = 0.4053g 
p(B4C) = —-0:8869g/em* 
N = pNA/ M «Eq. 2.59» 
N(B,C) = ———————=9.6679 10° 10” atoms/cm"-sec 
o( B4C) = 40,(B) + в((С) 
= 4(6,(B) + с;(В)) + (c4(C) + 6,(C)) 
= 4(759 + 3.6) + (0.0034 + 4.75) <Table II.3> 
= 3055.15b 
X(B4C)-N(B,C) o(B4C)-29.6679 10° 10^ 3055.15 10%” = 29.5369 cm” 
M(AI) = 26.9815 


N(AD = —————=6-0239 10? 107 atoms/cm"-sec 

o, ( AD) = o, (Al) + o,( AD) = 0.23 + 1.49 = 1.72b 

X(AD-2N(AD (АІ) = 6.0239 10? 107 1.72 10” = 0.1036 cm" 
N(Compressed Al) = ——————=5.8208 10? 10^ atoms/cm"-sec 


L (Compressed Al) = N(Compressed АІ) 6,(Al) = 5.8208 107 107 1.72 10^ 
= 0.1001 em” 
Penetrating probability 


-1 10° 


9. What is the probability that a neutron can move one mean free path without interacting in a medium? 


[sol] the probability that a neutron can move through this distance without having a collision 


= <page 58> 


= <опе mean free path, X, 21/2» 
= e! = 0.3679 


10. A wide beam of neutrons of intensity Оо is incident on a thick target consisting of material for which o, » Сз. 
The target area is A and its thickness is X. Derive an expression for the rate at which neutrons are absorbed in 
this target. 


[sol] neutrons absorption rate = ос МАХ (I : intensity of the beam, А: area, X : thickness) «Eq. 3.2» 


a wide beam of neutrons of intensity : Oy 

target area : A 

target thickness : X 

Ot = Oa «ба,» 0, 

neutrons absorption rate = 6,09NAX = O@pX,AX 


11. Stainless steel, type 304 having a density of 7.86 g/cm’, has been used in some reactors. The nominal 
composition by weight of this material is as follows: carbon, 0.08%; chromium, 19%; nickel, 10%; iron, the 
remainder. Calculate the macroscopic absorption cross-section of SS-304 at 0.0253 eV. 


[sol] X, = No, 


-2 (М,ба0) 
М; = (wipNa) / (100М;) «Eq. 2.63» 
М, ьо = 12.01115 «Table П.3» 
М ьо = 51.996 


nickel — 58.71 


Moron = 55.847 


Nearbon = (0.08 7.86 0.6022 1075/(100 12.01115) = 3.1526 10“ 107 atoms/cm? 
Nehromium = (19 7.86 0.6022 1075/(100 51.996) 21.7296 10? 10” atoms/cm? 
Мика = (10 7.86 0.6022 102)/(100 58.71) = 8.0622 10° 10” atoms/cm? 
Niron = (70.92 7.86 0.6022 102)/(100 55.847) = 6.0108 10? 107 atoms/cm? 
Oa(carbon) = 0.0034 b «Table П.3> 
Oa(chromium) = 3.19 
Oa(nickel) = 4.43 b 
Oa(iron) = 5. 
у, = (3.1526 10“ 0.0034) + (1.7296 107 3.1) + (8.0622 10° 4.43) 

+ (6.0108 107 2.55) = 0.2426 ст"! 


12. Calculate at 0.0253 eV the macroscopic absorption cross-section of uranium dioxide (U5O), in which the 
uranium has been enriched to 3 w/o in ^?U. The density of UO; is approximately 10.5 g/cm’. 


[sol] 33, -Х с, 

Atomic weight of U : 1/M, = 1/100(3/235.0439 + 97/238.0508) «Eq. 2.65» 

> M, = 237.9595 

Molecular weight of UO, = 237.9595 +2 15.9994 = 269.9583 

N = pN4/ M XEq. 2.59» 

М(ОО)) = (10.5 0.6022 1074) / 269.9583 = 2.3423 107? 

3% is not a volume percent > volume percent of U = 237.9595/235.0439 0.03 = 0.030372 

volume percent of U = 1 - 0.030372 = 0.969628 

У,(00,) N(UO; 5G,(UO)) 
= (2.3423 10”) {0.030372 680.8 + 0.969628 2.7 +2 0.00027) 1074} 
= 0.5457 cm! 


13. The compositions of nuclear reactors are often stated in volume fractions, that is, the fractions of Фе volume 


of some region that are composed of particular materials. Show that the macroscopic absorption cross-section 
for the equivalent homogeneous mixture of materials is given by 


Where and are, respectively, the volume fraction and macroscopic absorption cross-section of the ith 
constituent at its normal density 


[Sol] 


"where =N, =N, =N 


14. Calculate and (at 0.0253eV) for the fuel pellets described in Problem 2.63. The pellet density is 
about 10.6 g/c 


[sol] 
At Prob. 2.63 pellet composed of mixed oxides, U апа Pu with the Pu comprising 
approximately 30w/o of the mixture. 100% апа plutonium contains the following isotopes 


(70.5w/o), (21.3w/o), (5.5w/o) and (2.75w/o) 


—J=0-7- , )=0.3 


NO = 


NO = 


a( ) + а( ) + а( + а( 


0.7038 
=850.864840 


а( 
а(= 0.299 


7.056 =4.0818 


15. Using the fact that the scattering cross-section of 15 approximately 9 b from 0.01 eV to 200 eV, 


estimate the radius of the nucleus and compare with Eq.(2.3) 
[Sol] 

1) Atlow energy .. assume that the cross-section is 4 

50,4 =9b=9 


К- — — — 298.462843753 


2) using Eq 2.3 
R-1.25fm - Eq(2.3) 
В = 1.25 =7.418090175 cm7.42 


16. Using the Breit- Wigner formula. Compute and plot in the vicinity of the first resonance in ‚ Which 
occurs at an energy of 6.67eV. The parameters of this resonance are: =1.52meV(meV = millielectron volts), 

and g = 1. [Note: For reasons given in Section 7.3(see in particular Fig.7.12), the value of 
computed from the Breit- Wigner formula do not concide with those measured with targets at room temperature. | 


[Sol] 
Breit- Wigner Formula 


g=1, =1.52meV, , = 1.52 + 26 =27.52meV 


2.039x 








-0.00 





17. Demonstrate using the Breit- Wigner formula that the width of a resonance at half its height is equal to . 


[Sol] 


— (Ё-)-0 So, 


18. There are no resonances in the total cross-section of from 0.01 eV to cover 1 Me V. If the radiative 
capture cross-section of this nuclide at 0.0253 eV is 3.4 mb, what is the value of at 1 eV? 


[sol] , it can be written as 


SO 


19. The first resonances in the cross-section of aluminum, which is due entirely to scattering, occurs at 5.8keV. 
The absorption cross-section at 0.0253 eV is 0.23b. Calculate for 100 eV: 


(a) 
[Sol] 


—, SO 


—=0:23b ——=3.658373956 


(b) 

[Sol] 

Low energy so it can assume =4=4 

So = =1.767145868 
(c) 


[Sol] 


20. Calculate : 
(a) water of unit density at 0.0253 eV; 


[Sol] 
А at 0.0253eV 


N =— 


(b) water of density 0.7g/ at 0.0253eV; 


[Sol] 


ш--шт — —_234— 


(c) water of density 0.7g/ at 1 eV. 


[Sol] 


2]. The first resonance in the scattering cross-section of the nuclide occurs at 1.24 MeV. The separation 
energies of nuclides , , апа аге 7.00, 7.50, and 8.00 Me V, respectively. Which nucleus and at 
what energy above the ground state is the level that gives rise to this resonance? 


[Sol] 


9 9 


. The first resonance in the scattering cross-section of the nuclide occurs at 1.24 MeV 


So it should occur the the energy is 
8.00MeV + 1.24MeV = 9.24MeV 


22. There is a prominent resonance in the total cross-section of аг 646.4keV. At what energy, measured 
from the ground state, is the energy level in that corresponds to this resonance?[Hint: Use the masses of 
neutral and to compute the binding energy of the last neutron in ] 


[Sol] 


The energy level is 7.646MeV + 646.4 =8.2964MeV 


23. The excited states of | occur at the following energies (in Me V) measured from the ground state : 0.871, 
3.06, 4.55,5.08,5.38,5.70,5.94, etc. At roughly what energies would resonances be expected to appear in the 
neutron cross-section of 


[Sol] 


The first energy above the 4.144MeV is 4.55 MeV s 
4.55 — 4.144 -0.406МеУ this energy will be appear when the resonances 


24. Using Eq. (3.28), compute and plot E'/E as a function of angle from 0 to л for A= 1, 12 and 238. 


[Sol] Eq. (3.28) —— 
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25. A 2-MeV neutron traveling in water has a head-on collision with an '°O nucleus. 


(a) What are the energies of the neutron and nucleus after the collision? 


[Sol] ——;Ж=2 MeV, A-16, с=л 


(b) Would you expect the water molecule involved in the collision to remain intact after the event? 


[Sol] By the laws of conservation of momentum, the water molecule was moved towards direction of 
initial neutron. 


26. А 1-MeV neutron strikes а "C nucleus initially at rest. If the neutron is elastic scattered through an angle of 
90°: 


(a) What is the energy of the scattered neutron? 


[Sol] ——-E-1 MeV, A=12, в=л/2 


(b) What is the energy of the recoiling nucleus? 
[Sol] 
(c) At what angle does the recoiling nucleus appear? 


[Sol] 


Y-axle) 
sin 


Other solution 





X-axle) p 
cos 


27. Compute and plot the average fractional energy loss in elastic scattering as a function of the mass number of 
the target nucleus. 


[Sol] Eq. (3.32) — - —— 
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28. Show what the average fractional energy loss in% in elastic scattering for large А is given approximately by 


[Sol] — - — —— — - — 


29. А 1.5-Ме V neutron in a heavy water reactor collides with an ^H nucleus. Calculate the maximum and 
average changes in lethargy in such a collision. 


[Sol] Eq. (3.28) — 


At б-л, Lethargy is maximum. 
A=2, — - 


Average change, — — - 


30. Suppose that a fission neutron, emitted with energy of 2 MeV, shows down to an energy of 1 eV as the result 
of successive collisions in a moderator. If, on the average, the neutron gains in lethargy the amount & in each 
collision, how many collisions are required if the moderator is 


(a) Hydrogen, 


[Sol] - 
After the number of n time collisions was occurred, - 


(b) graphite? 


[Sol] - 


After the number of n time collisions was occurred, = 


31. The 2,200 meters-per second flux in an ordinary water reactor is 1.5 х 10" neutrons/cm7-sec. Ate what rate 
are the thermal neutrons absorbed by the water? 


[Sol] 


32. At one point in a reactor, the density of thermal neutrons is 1.5 x 10? neutrons/cm’. If the temperature is 
450°C, what is the 2,200 meters-per-second flux? 


[Sol] 


33. A tiny beryllium target located at the center of a three-dimensional Cartesian coordinate system is 
bombarded by six beams of 0.0253-eV neutrons of intensity 3 x 10° neutrons/cm’-sec, each incident along а 
different axis. 

(a) What 1s the 2,200 metersOper-second flux at the target? 

[Sol] 

(b) How many neutrons are absorbed in the target per cm^-sec? 


[Sol] 


34. When thermal neutrons interact with "N, what is the probability that absorption leads to radioactive capture? 
[Sol] — — 


35. The control rods for a certain reactor are made of an alloy of cadmium (5 w/o), indium (15 w/o), and silver 
(80 w/o). Calculate the rate at which thermal neutrons are absorbed per gram of this material at a temperature of 
400°C in a 2,200 meters-per-second flux of 5510" neutrons/cm.-sec. [Note: Silver is a 1/v absorber.] 


[sol] Eq.(3.42)- (3.44) 


To calculate the atomic abundance of each elements, 








36. From the data in Table 3.3, would you expect ^" 


expect fission to be possible? 


Th to be fissile? If not, at what neutron energy would you 


[sol] ТЬ + їп = ГВ, the critical energy of ^?Th is 5.1Me У, and B.E is 6.5Me V. 
К.Е = critical energy - B.E = 1.4Me V 


г. Higher than 1.4Ме V 


37. Two hypothetical nuclei, ^7. and “*'Z, of atomic weights M(^Z)2241.0600 and M(^*'Z)2242.0621 have 
critical fission energies of 5.5 MeV and 6.5 MeV, respectively. Is the nuclei ^Z fissile? 


[sol] 


38. Fission can be induced when y-rays are absorbed by a heavy nucleus. What energy y-rays are necessary to 
induce fission in 

[sol] 5.75 Me V 

[sol] 5.85Me V 

(c 9Pu 

[sol] 5.5Me V 

39. Cross-sections of ^?U at 1 MeV are as follows: 0,74.0b, б::1.40, o=1.2b, о.=1.3Ъ. The cross-sections for 
neutron-producing and charged particle reactions are all negligible. Compute at this energy 


(a) the total cross-section 


[sol] SO, 


(b) the capture-to-fission ratio а 
[sol] - — — 


40. The fission cross-product BIT has a half-life of 8.05 days and is produced in fission with a yield of 2.9%-that 
is, 0.029 atoms of Ч are produced per fission. Calculate the equilibrium activity of this radionuclide in a 
reactor operating at 3,300 MW. 


[sol] — 


41. Fission-product activity measured by at the time t; following the burst of a nucleus weapon is found to be бо. 
Show that the activity at the time t = 7"to is given approximately by a= 0/10". This is known as the 7-10 rule in 
civil defense. 


[sol] At 83p, fission product activity 


42. Suppose that radioactive fallout from a nuclear burst arrives in a locality a hour after detonation. Use the 
result of Problem 3.41 to estimate the activity 2 week later. 


[sol] 


43. The yields of nuclear weapons are measured in kilotons (KT), where 1 KT —2.6*10? MeV. With this in mind, 
(a) How much 2350 is fissioned when a 100-КТ bomb is expoded? 


[sol] 


(b) What 1s the total fission-product activity due to this bomb 1 min, 1 hr, and 1 day after detonation? [Note: 
Assume a thermal energy release of 200 MeV per fission.] 
[sol] 

for 1 min, == pone 


for 1 hr, -- 
for 1 day, 


44. A research reactor is operated at a power of 250 kilowatts 8 hours a day, 5 days a week, for 2 years. A fuel 
element, 1 of 24 in the reactor, is then removed for examination. Compute and plot the activity of the fuel 
element as a function of time up to 2 years after removal. 


[sol] Fission product activity — 











300000 - 
200000 
activity(Ci) 
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45. The spontaneous fission rate of ^*U is 1 fission per gram per 100 sec. Show that this is equivalent to a half- 
life for fission of 5.5*10P years. 


[sol] ———— 


46. Compute and plot the parameter 1) at 0.0253 eV for uranium enriched in ^?U as a function of its enrichment 
in weight percent ^?U. 


[Sol] 


N; = (х/100) - (pni /M;) where x = w/o, p 719.1 g/cm? 


== у(235)• /[ + ] 
= ү(235)•3(235)•с(235) / N(235)*o, (235) + N(238)*o, (238) 


0.6889X/ [0.333X+(100-X)1.305 x 107] 


= 0.6889X/ (0.1305+0.332Х) = 2.075(Х+0.332-0.332)/(Х+0.393) 
= 2.075 — 0.6889/(X+0.393) 


47. Suppose that 1 kg of U undergoes fission by thermal neutron. Compute the masses (or mass equivalents) 
in grams for the following, which are produced: 

(a) neutrons 

(b) B-rays 

(с) y-rays 

(d) neutrinos, 

(e) kinetic energy 

(f) fission products. 


[Sol] 
A number of fission occurred by 1kg of ^^U 


1000 x 0.6022 х 1074/(235.0439) = 2.56 x 107^ fissions 


(a) neutron 
v(235)=2.418 


2.418 x 2.56 x 10” neutrons 


(b) B-rays 
2.56 x 10 В-гауѕ x (8 Me У/р-гау) = 2.408 x 107 Me V 


(c) y-rays 
2.56 x 10 x 72 1.792 x 10^ Me V 


d) neutrinos 
2.56 x 105 х 12 = 3.072 x 10? Me V 


(e) kinetic energy 
2.56 x 10% x 5 =1.28 x10" Me V 


(f) fission products. 


2.56 x 10“ x (84-7412) = 6.912 x 10° Me V 


48. The reactor on the nuclear ship Savannah operated at a power of 69 MW 
(a) How much U was consumed on a 10.000-nautical mile voyage at an average speed of 20 knots? 
(b) This is equivalent to how many barrels of 6.5-millon-Btu/barrel bunker-C oil? 


[Sol] 

a) Consumption rate = 1.05(1+а)-Р (g/day) 
=1.05(1+0.169)-69 a value can be seen from Table 3.4 
=84.694 (g/day) 


Hours under voyage 
t = 10000 nautical mile/20knots = 500Hr = 20.83days 
Masses of consumed 20 


= 84.694 х 20.83 = 1764.46 © 


b) Total work = 69MW х500Нг 
= 34500 MW-Hr x (1Btu/2.931 x 10" MW-Hr) 
= 1.777 x 10! Btu 


The amounts of consumed bunker-C oil 
21.777 x 10" Btu / (6.5 x 10° Btu/barrel) 21.811 x 10* barrels 


49. Consolidated Edison's Indian Point No 2 reactor is designed to operate at a power of 2.758 MW. Assuming 
that all fission occur in ^"?U, calculate in grams per day the rate at which ^?U is 

(a) fissioned, 

(b) consumed. 


[Sol] 
a) Fission rate 


= 2758 MW х (10° Ј/ MW) x (fission/200 Me V) x (Me V/1.6 x 10? J) x (86400 sec/day) 
—7.447 x 10" fissions/day 
b) Consumption rate 
= 7.447 х 10" fissions/day x 235 g/mole 0.6022 x 10" fissions/mole x (140.169) 
= 3397.2 g/day 


50. Referring to the preceeding problem, what is the total accumulated activity of the fission products in the 
Indian Point No. 2 reactor 1 day after shutdown following 1 year of operation? 


[Sol] 
a 21.4 x 105 x P [t? — (t T) ?] 


-14 x10°x 2758 x [1?— (14365) ?] 
= 2.675 x 10 С1 


51. The photoelectric cross-section of lead at 0.6 Me V is approximately 186. Estimate бүс at this energy for 
Uranium. 


[Sol] 
Ope ~Z" 
Сре = const. Z^ 

From fig. 3.14, the value of n for 0.6 Me V-y-ray 4.35 
18b-const.(82)*^ 


Const. = 18 x (82) = 8.515 x 10°C 
For Uranium 
Gpe = 8.515 х 10?Z^ 


= 8.515 х 10°(92)*** 
= 29.69b 


52. A 2-Me V photon is Compton scattered through an angle of 30°. 
(a) What 18 its energy after scattering? 

(b) What is the recoil energy of struck electron? 

(c) At what angle does the electron appear? 


[Sol] 
а) Е' = E E. /2(1-COSo) *E,2(2x0.511) [2х (1-COS30°)+0.511]=1.312 Me V 


b) Recoil energy = E- E'=2-1.312=0.688 Me V 


= о 
c) 6-30 
E. p' ний” 
227 Scattered neutron 
Г, р = ; 9 
O [7 SAN СЕ Кы Е 
iS _ x 
Incident [ир ^-: 9 
neutron nucleus м 


а x Kecoiling nucleus 
E.P нар 
Psing=p'sin30° 
ѕіпф=р'/Р x 1/2 


p'=h/4' = h/[(1.240 x 109) Е!= 4.135710 eV-sec/[(1.240 x 10^) 1.312 x 104]-cm 
-4.376 х 10" Me У:ѕес/ст 

= (4.376 х 107 Me У:ѕес/ст) х [(1.6 x 10* erg)/Me V] x [(g-cm?/sec?y/erg] 

= 7.01 x 1077 g-cm/sec 


Р=т,-Е,'= V2 x 9.10956 x 107 x 0.688 x 1.602 x 1024.48 x 1077 g-cm/sec 


sing= p'/2P=7.01 x 1077/4.48 x 1077-0.780 
ф=51.44° 


53. Show that eq. (3.62) follows from eq. (3.61). 


[Sol] 
Е' =EE,/2(1-COSo) +E, 
À = h/p = hc/E , E-hc/ À 
he/ 4' = (he/ А. - һс/ А) /[hc/ А (1-COSo)+ hc/ 2,)] 
he/ A'[he/ X (1-COSo)* he/ А,)]=( Һе)2/2 Ae 
(1-COSo) /А'-А + (1/^/ №) =1/K * А) 
1/A.-(1/A-1/2') = (1-COSo) /.-. 
1/) = 1/( A'- А (1-COSo) 
À'- A= 2,(1-COSo) 


54. Show that the minimum energy of the scattered photon in Compton scattering is given by 


(Е) min = ЕЕ. /2Е +E, 
And that for E» » E,, 


(E) „= E,/22 0.255 Me V. 

[Sol] 
a) E'=EE,/2(1-COSo) +E, 

At o=7 E' = (E!) min 

E' = EE'/2(1-COS л) +E = EE'/2E +E, 
b) For EE, 

2Е-Е, 2E 
(E) mn = EE, /2E=E,/2 = 0.255 Me V 


55. What is the minimum energy of a Compton-scattered photon if its original energy is 
(a) 0.1Me V, 
(b) 1 Me V, 
(c) 10 Ме V? 


[Sol] 
a) 0.1 Me V 


(E')min= = 0.1 x 0.511) x [2 x 0.1)+0.511] = 0.072 Me V 
b) 1 Me V 

(E')min=(1 х 0.511) x [(2 x 0.1)20.511]- 0.204 Me V 
c) 10 Me V 

(ЕЭ) nin=(10 x 0.511) x [(2 x 0.1)+0.511]= 0.249 Me V 


56. Calculate the mass attenuation coefficient of silica glass ($105, р = 2.21 g/cm”) for 3- Me V y-rays. 


[Sol] 
М(51О:) = 28.089+(2 x 15.9994) = 60.0848 


(u/p)s; = 0.0367 cm?/g 
(u/p)o= 0.0359 cm?/g  ,for3 Me V- photon 


u/p = Y';co;(u/ р) = (28.086/60.0848)-0.0367 --[ (2 х 15.9994)/60.0848]-0.0359 
= 0.0363cm"/g 


57. derive Eq. (3.69) 


[sol] 


58. The mass attenuation coefficient of lead at 0.15MeV is 1.84сш /в. At this energy, the principal mode of 


interaction is by the photoelectric effect. What thickness of lead is required to reduce the intensity of a 0.15- 
MeV y-ray beam by a factor of 1000? 


[sol] Table II.3 


59. The density of air at standard temperature and pressure (0°C and latm)is g/cm^. Compute the 
man free paths of photons in air under these conditions and compare with the corresponding mean free paths in 


unit-density water at the following energies: 
(a) 0.1 MeV 


[sol] (1) Air 
Table II.4 - 


(2) Water 
Table II.4 - 


(b) 1 MeV 


[sol] (1) Air 
Table II.4 - 


(2) Water 
Table II.4 - 


(c) 10 MeV 
[sol] (1) Air 
Table II.4 - 


(2) Water 
Table II.4 - 


60. At 1 MeV, the Compton cross-section per electron is 0.2112 b, and the Compton energy absorption cross- 
section per electron is 0.0929 b. 
(a) What 1s the average energy of the recoiling electron in a Compton interaction at this energy? 


[sol] 


(b) Compute the Compton mass attenuation and mass absorption coefficients at 1 MeV for (1) aluminum, (11) 
water. 


[sol] (1) Aluminum 
Table II.3 М-26.9815 


(2) Water 
Table II.3 M=18.0153 


61. A beam of 0.1-MeV y-rays with an intensity of  y-rays/cm'-sec is incident on thin foils of (i) 
aluminum, (ii) water. At this energy, the Compton cross-section per electron is 0.4929 b, and the Compton 
energy absorption cross-section per electron is 0.0685 b. Calculate the energy extracted from the beam per unit 
volume of the foils due to 

(a) Compton scattering 


[sol] (1) Aluminum 
Table II.3 M=26.9815, 


(2) Water 
Table II.3 M=18.0153, 


(b) the photoelectric effect 


[sol] (1) Aluminum 
Assume : 


? 


— — — —20-04703 


(2) Water 


— — —— 270-0024 


62. The absorption of radiation 18 often measured in units called rads, where 1 rad is equal to the absorption of 
100 ergs per gram. What intensity of IMeV y-rays incident on a thin slab of water is required to give an 
absorption rate of 1 rad per second? 


[sol] 1rad-100erg/g- 
(Absorption rate)- =EI — 


63. Determine the range of 5-MeV a-particles in the following media: 
(a) air at 15°С, latm; 


[sol] Figure 3.23 5MeV 


(b) aluminum; 
[sol] Table ПЗ Mz26.9815 


(c) lead; 


[sol] Table II.3 M-207 .19 


(d) unit-density water; 


[sol] Table II.3 M=18.0153 


(e) air at 300°C, 10atm; 


[sol] (0: Air at 15°C, latm, © :Аіг at 300°C, 10atm 





64. Determine the relative stopping powers of the media in the preceding problem. 


[sol] (a) Air 


(b) Aluminum 


(c) Lead 


(d) Unit-density water 





(e) Air at 300°C, 10atm 





65. Compare the apparent mass attenuation coefficient of 2-MeV (maximum energy) В-гауѕ with the mass 
attenuation coefficient of 2-MeV y-rays in aluminum. 


[sol] - — 


66. Compare the maximum ranges of 3-MeV a-rays and p-rays in air at standard thmperature and pressure. 


[sol] 





(At 20°С, 101.325kPa) 


67. Near the surface of a flat fuel element in an operating reactor, fissions are occurring at the constant rate of S 
fissions/cm"-sec. Given that the average range of the fission fragments is В, show that the rate at which such 
fragments would escape рег сш /вёс from the surface of the fuel if it were not clad is equal to SR/2. 


[sol] 


R 


(+) 
(№) 





Flat fuel 


(The rate at which such fragments would escape per ) 


-SRfissions =@+@ 


Chapter 4 


Nuclear Reactors and 
Nuclear Power 


1. If 57.596 of the fission neutrons escape from a bare sphere of ^?U, what is the multiplication factor of the 
sphere? In this system, the average value ofn is 2.31. 


[sol] ———— 


2. Measurements on an experimental thermal reactor show that, for every 100 neutrons emitted in fission, 10 
escape while slowing down and 15 escape after having slowed down to thermal energies. No neutrons are 
absorbed within the reactor while slowing down. Of those neutrons absorbed at thermal energies, 6090 are 
absorbed in fission material. (a) What is the multiplication factor of the reactor at the time these observations are 
made? (b) Suppose the thermal leakage is reduced by one third. How would this change the value of А? [Note: 
The values of 7 and v for the reactor fuel are 2.07 and 2.42, respectively.] 


[sol] (a) 
(b) К 2— ——— 


3. (a) Show that the energy released in the nth generation of a fission chain reaction, initiated by one fission, is 
given by 


? 


where К is the multiplication factor and Eg is the recoverable energy per fission. (b) Show that the total energy 
released up to and including the nth generation is given by 


[sol] (a) One generation : 
Two generation : 


'. nth generation : 


(b) Total energy : 


4. Show that the fraction, F , of the energy released from a supercritical chain reaction that originates in the final 
m generations of the chain is given approximately by 


? 


provided the total number of generations is large. 


[sol] 


5. (a) Most of the energy from a nuclear explosion is released during the final moments of the detonation. Using 
the result of the previous problem, compute the number of fission generations required to release 99% of the 
total explosive yield. Use the nominal value  . (b) If the mean time between generations is the order of 

sec over what period of time is energy released during a nuclear explosion? 


[sol] (a) — = = 
0.99= - — | —k-2 


9 


(b) (Period of time)= 


б. A burst of neutrons from a pulsed accelerator is introduced into a subcritical assembly consisting of 

an array of natural uranium rods in water. The system has a multiplication factor of 0.968. Approximately 80% 
of the incident neutrons are absorbed in uranium. (a) How many first-generation fissions do the neutrons 
produce in the assembly? (b) What is the total fission energy in joules released in the assembly by the neutron 
burst? 


[sol] (a) 


(b) 
—=6.01917J 


7. It is found that, in a certain thermal reactor, fueled with partially enriched uranium, 13% of the fission 
neutrons are absorbed in resonances of О and 3% leak out of the reactor, both while these neutrons are 
slowing down; 5% of the neutrons tat slow down in the reactor subsequently leak out; of those slow neutrons 
that do not leak out, 82% are absorbed in fuel, 74% of these іп 7°U. (a) What is the multiplication factor of this 


reactor? (b) What is its conversion ratio? 


[sol] 
(a) 


(b) 


9. A natural uranium-fueled converter operates at a power of 250MWt with a conversion ratio of 0.88. At what 
rate is “Ри being produced in this reactor in kg/year? 


[sol] А 


— Table 3.4 


73.59kg/year 
9. Assuming a recoverable energy per fission of 300 Me V, calculate the fuel burnup and consumption rates in 
g/MWd for 
(a) thermal reactors fueled with 233U or 239Pu 


[sol] (1) case of U-233 


Fission rate = 2.88 * 10?'P  fissions/day (replace 200 Me V to 300 Me V) 


—1+ = = 


Consumption гаї : 1.1382(1.087)=1.1355 g/MWd 


(2) case of Pu-239 


Burn-up rate : — 
Consumption rate : 1.142(1-0)Р g/day 


1+0=— — — 50 consumption rate is 1.4593 g/day 


(b) fast reactors fueled with 239Pu. | Хос: In part (b), take the capture to fission ratio to be 0.065.] 


[sol] Burn up rate : — 


Consumption rate : 
Devide by unit power[M W] 


10. Because of an error in design, a thermal reactor that was supposed to breed on the 232Th-233U cycle 
unfortunately has a breeding ratio of only 0.96. If the reactor operates at a thermal power level of 500 megawatts, 


how much 232Th does it convert in 1 year? 
[sol] in 120page 


Consumption rate 


C(breeding ratio)=0.96, in 1 yr 


11. What value of the breeding gain is necessary for a fast breeder operating on the 238U-239Pu cycle to have 
an exponential doubling time of 10 years if the specific power for this type of reactor 18 0.6 megawatts per 


kilogram of 239Pu? 


[sol] Consumption rate per unit power [per MW] 


12. Per GWe-year, a typical LMFBR produces 558 kg and consumes 789 kg of fissile plutonium in the core 
while it produces 455 kg and consumes 34 kg of fissile plutonium in the blanket. What is the breeding ratio for 
this reactor? [Note: As would be expected, there is a net consumption of plutonium in the core and a net 
production in the blanket, with a positive output from the reactor as as whole. By adjusting the properties of the 
blanket, it is easy to make breeders into net consumers of plutonium if dangerously large stockpiles of this 
materials should ever accumulate in the world. | 


[sol] —————————5 


13. A certain fossile fueled generating station operates at a power of 1,000 MWe-at an overall efficiency of 38% 
and an average capacity factor of 0.70. 


(a) How many tons of 13,000 Btu per pound of coal does the plant consume in 1 year? 


[sol] ——_— Е 


(b) It an average coal-carrying railroad car carries 100 tons of coal, how many саг loads much be delivered to 
the plant on an average day? 


[sol] — —— 


(c) If the coal contains 1.596 by weight sulfur and in the combustion process this all goes up the stack as SO2. 
how much SO2 does the plant produce in 1 year? 


[sol] 


14. Had the plant described in Problem 4.13 been fueled with 6.5 million Btu per barrel bunker-C fuel oil 
containing .3790 sulfur, 


(a) how many barrels and tons of oil would the plant consume in 1 year ? 
[sol] 1 barrel=117.35 liter 

Density of Bunker C oil = 0.96 g/ml 

112.7 kg/barrel 


(b) how much SO2 would the plant release in 1 year? [Note: 1 U.S. petroleum barrel=5.61 cubic feet=42 U.S 
gallons; the density of bunker-C oil is approximately the same as water. | 


[sol] 


15. Rochester Gas and Electric’s Robert Emmett Ginna nuclear power plant operates at a net electric output of 
470 megawatts. The overall efficiency of the plant is 32.3%. Appoximately 60% of the plant’s power comes 
from fissions in 235U, the remainder from fissions in converted plutonium, mostly 239Pu. If the plant were 


operated at full power for 1 year, how many kiomgrams of 235U and 239Pu would be 
(a) fissioned? 


[sol] 
Burn up rate of each 


2) Ph-239 


(b) consumed? 


[sol] 1) U-235 


2)Ph-239 


16. Consider two electrical generating stations-one a fossil fuel plant and the other nuclear-both producing the 
same electrical power. The efficiency of the fossil fuel plant is (eff), whereas that of the nuclear plant is (eff),. 
(a) Show that the ratio of the heat rejected to the environment by the two plants is given by 


[sol] Qr : thermal energy output from the heat source. 
Ос : heat energy rejected to the environment 


In the exercise, W = (eff; х Ор; = (eff, x Оһ, 


(b) Evaluate О,/О,, for the case where (eff); = 38% and (eff), = 33% 
[sol] — ---- --- — — 


17.Ап MSBR power plant produces 1,000 MWe at an overall efficiency of 40%. The breeding ratio for the 
reactor is 1.06, and the specific power is 2.5 MWt per kilogram of ^"U. 
(a) Calculate the linear and exponential doubling times for this reactor. 


[sol] 


(b) What is the net production rate of ^?U in kg/year? 


[sol] шин — 


18. A 3,000MW reactor perates for 1 year. How much does the mass of the fuel change during this time as the 
result of the energy release? 


[sol] The fissioning of 1.05g of U yields IMWd. 


19. Referring to the nominal LWR cycle described in Fig. 4.37, The plant operates at an overall efficiency of 


33.4%,compute 
(a)The specific burnup of the fuel in MWd/t 


[sol] — 


(b) The fractional burnup of the fuel 

[sol] —— 

(c) The enrichment of Фе fresh fuel 

[sol] —— 

(d) The enrichment (235U) of the spent fuel 
[sol] —— 


(e) The fraction of the power originating in fissions in ^?U and plutonium, respectively 


[sol] fissioning of 1.05g of 20 yields IMWd. 


fissioning of 1.052 of “Ри yields IMWd. 


(f) The 30-year requirement for uranium for this system. 


[sol] 150,047kg of 0308 is converted to UF6 in year. 


22. Using the data in Fig. 4.41 [Note : The reactor operates at an efficiency of 30%.], compute 
(a) The enrichment of the fresh fuel 


[sol] 


(b) The residual enrichment of the spent fuel 


[sol] mass of 2550) = 
A mount of 60% of neutrons are absorbed in fission material. 


(c) The fraction of power originating in fissions in "0 


[sol] — 
To yield 2500M We-yr, —— —— мы 





Fraction= 
(d) The burnup of the fuel. 
[sol] АЕ о үн МИН 


23. Reliable nuclear weapons cannot be made using plutonium containing much more than 7 w/o Ри because 
this isotope has a high spontaneous fission rate that tends to pre-initiate - that is, fizzle — the device. Using Fig. 
4.37 and Fig. 4.39 through4.40, determine whether any of these commercial power systems produce weapons- 


grade plutonium (<7 w/o Ри). Assume that all non-fissioning plutonium 18“. Ри (Some of it is actually ^" Pu). 


[sol] 


24. Compute and plot annual and cumulative uranium requirements through the year 2015 for the projected 
nuclear electric capacities given in Table 4.8 for each of following assumptions: (a) all reactors are PWRs with 
no recycling; (b) all reactors are PWRs, and all operate with recycling after 2000; (c) all reactors are PWRs until 
2000 after which 20% of new reactors are LMFBRs. [note: for simplicity, ignore the extra fuel required for 
startup; that is, take the annual uranium needs to be the values in Table 4.7 divided by 30.] 


[sol] According to Table 4.7 Natural U requirements for 1 yr in LWR 
Once through: 4260/30- 142t 
Recycling: 2665/30=90t 

For once through cycle, 


a) 





| 6.49 | 
2000 10.11 


2005 14.02 
2010 18.21 
2015 22.54 


Breeding ratio(cumulation) 


10.11 


12.55 


15.17 
17.88 


c) In LMFBR, 36/30=1.2t/yr for U-Pu recycling 


Breeding ratio Breeding ratio(cumulation) 
1985 





1990 
195 |  — 178 | — 406  «». — 





2000 





25. A nuclear fuel-fabricating company needs 10,000 kg of 3 w/o uranium and furnished natural uranium (as 
ОЕ6) as feed. Assuming a tails assay of 0.2 w/o, (a) how much feed is required? (Б) How much more will the 
enrichment cost? 


[sol] 
a) 


Product 






M p^ X p 







Enrichment 
plant Tails 


XT 


Feed 
M p Xp 


Хр= 0.00711, Хт= 0.002, Xp- 0.03 
When the 10,000 kg of U is M,, mass of UO; is 


10,000 x [(238+2 x 16)/238]=11345 kg 


My = M, +Mp 
ХеМе=ХтМт+ХрМр 
— My -z(Xp-X4)/ (Xr-Xr) x Мь-54795 kg 
=(0.03-0.002)/(0.0071-0.002) х 10000 = 54795 kg 


In ОЕ, 54795 x [(238+19 x 6)/2381-81041 kg 

b) 

y(0.03)2(1-2-0.03)1n[(1-0.03)/0.03]23.267 5 

y(0.00711)= (1-2-0.00711)In[(1-0.00711)/0.00711]-4.869 
ү(0.002)= (1-2-0.002)In[(1-0.002)/0.002]=6.188 


SWU= Mp -y(0.03)+ M, -y(0.002)- Ме -y(0.00711) 


= 10000-3.2675+(54795-10000)-6.188 - 54795-4.869 
= 43070 kg 


$130.75/SWU 
Enrichment cost= 130.75-43070=$5,631,400 


26. Suppose that the tails enrichment in the preceding problem were decreased to 0.15 w/o. (a) How much 
uranium feed would this save? (b) How much more would the enriched fuel cost? (c) What is the cost per kg of 
feed saved? 


[sol] 
a) 


M' = (ХЬ-Х'т)/ (Xg-X' 4) x Mp 
=(0.03-0.015)/(0.0071-0.015) х 10000 = 50802 kg 

My - M'p=54795-50802= 3993kg 

In ОЕ, 3993 . [(238+19 х 6)/238]- 5906kg of ОЕ, ore saved 


b) 
y(0.03)2(1-2-0.03)1n[(1-0.03)/0.03]23.267 5 


y(0.00711)= (1-2-0.00711)In[(1-0.00711)/0.00711]-4.869 
ү(0.015)= (1-2-0.015)In[(1-0.015)/0.015]26.4813 


SWU= Mp -y(0.03)+ Мт -y(0.015)+ Mr -y(0.00711) 


= 10000-3.2675+(54795-10000)-6.4813 - 54795-4.869 
= 49770kg 
Enrichment cost = 130.75-49770=$6,507,428 


$6,507,428-$5,631,400=$876,028 
c) 876,028/3,993=$219.39 


27. If a customer supplies DOE as part of his feed-uranium that is partially enriched above (or depleted below) 
the natural uranium level-the total amount of natural uranium feed that he must furnish to obtain a given 
amount of product is reduced (or increased) by an amount equal to the natural uranium feed needed to provide 
the partially enriched (or depleted) feed. The customer receives a credit (or debit) for the separative work 
represented by his enriched (or depleted) feed. 


Suppose that the fuel-fabricating company in problem 4.25 offers to supply as part of its feed 10,000 kg of 1% 
assay uranium. How much natural uranium feed is required in addition, and how much will the total job cost? 


[sol] 


a) 
The amount of Mp generated from the 100000f 1 w/o U 


Mp, = (Xr-Xr )/ (XPp-X7) х Ме 


=(0.01-0.002)/(0.0071-0.002) x 10000 = 2857 kg 
The amount of 3 w/o U produced by the Natural U 
10000-2857-:7143Ко 


Natural U needed to make the 7143kg of 3 w/o U 
М=(ХЬ-Хт )/ (Xy-X1) х Mp 


=(0.03-0.002 )/ (0.00711 -0.002 ) x 7143 239139.7kg 
b) 
ү(0.01)=4.503 
ү(0.00711)= 4.869 
ү(0.002)=6.188 
ү(0.03)= 3.267 


i) For 1 w/o О 
SWU= Mp -y(Xp)+ Му Xr) -Mr Y (Xp) 


= Mply(Xp)-y(X7)] - МЕҮ(Хр)- YCX1)] 
=2857(3.267-6.188)-10000(4.503-6.188) 
=85047 kg 
ii) For natural U 
SWU- 7143(3.267-6.188) -39139.7(4.869 -6.188) 
= 30760.6kg 


8504.7+30760.6=39265.3k¢g 


Enrichment cost=$130.75 x 39265.3=$1,876,881 


28. Repeat problem 4.27 for the case in which the customer furnishes 10,000 kg of 6% feed. 


[sol] 
a) 


The amount of Mp generated from the 10000 of 0.6 w/o U 
Мы = (Хк-Хт)/ (ХРь-Хт) x Ме 
=(0.006-0.002)/(0.03-0.002) x 10000 = 1429 kg 


The amount of 3 w/o U Produced by the Natural U 
10000-1249=8571kg 


Natural U needed to make the 7143kg of 3 w/o U 
М=(Хь-Хт )/ (Xy-X1) х Mp 


=(0.03-0.002 )/ (0.00711-0.002 ) х 8571 =46967kg 


b) 

ү(0.006)=5.049 
ү(0.00711)= 4.869 
ү(0.002)=6.188 
ү(0.03)= 3.267 


i) For 1 w/o О 
SWU- Mp "y (Xp) Мт : ү(Хт) -Mr (Xr) 
= Mp[y(Xp)-y(X7)] - Mely(Xy)- үСХт)] 
=1429(3.267-6.188)-10000(5.049-6.188) 
=7217 kg 
ii) For natural U 
SWU= 8571(3.267-6.188)-46967(4.869-6.188) 
= 36912 


72174369124-44129kg 


Enrichment cost=$130.75 x 44129=$5,769,866 


29. Referring to problem 4.28, suppose exactly 33,000 kg of 3.2 w/o UO2 15 required per reload. (a) How much 
UF6 must be given to the enrichment plant, assuming 1% loss in fabrication? (b) How much yellow cake must 
be used in the conversion to UF6 assuming 0.5% loss in conversion? 


[sol] 

a) 

Mass of U in UO, 

= 238/(238-16:2) :33000-29088.9Ко 


Mass of О in ОЕ, after enrichment, M, 
M,=1.00/0.99-29088.9=29382.7kg 


Mr = (Xp-Xq y (Xp-X4) X Mp -54795 Ко 
=(0.032-0.002)/(0.0071-0.002) х 29382.7kg = 172500kg 


Mass of ОЕ, given to enrichment plant 
=(238+19-6)/238 х 1725002255126kg 


b) 
U50; — Conversion to UF, (0.5% loss of total U) — ОЕ, 


Mass of U in 030; 
=(1.00/0.995) х 172500=173367kg 


Mass of О:О,М 


=[(3 -238+16 -8)/(3 -238)] -173367=204447kg 


30. It is proposed to produce 25 kg of 90 w/o for a nuclear weapon by enriching 20 w/o fuel from a research 
reactor. (a) How much fresh reactor fuel would be required? (b) Compute the total SWU required. (c) Compare 
the SWU/kg required to produce 20 w/o fuel starting from natural uranium with the SWU/kg for 90 w/o material 
beginning at 20 w/o. [Note: Assume the tails enrichment to be 0.2 w/o.] 


[sol] (a) My = ОМс-БЕа. 4.33» 
Xp = 0.9, хт = 0.002, xp = 0.2, Mp = 25 kg 
Мк = 0—25-kg = 113.384 kg 
(b) SWU = МЫ У (хр) — V(x7)] — Me[V(xp) – У(хт)] «Eq. 4.37» 
V(x) = (1-2х)  — <Eq. 4.35> 
V(xp) = V(0.9) 2(1— 1.8) —= 1.758 
V(x4) = V(0.002) = (1— 0.004)  —-6.188 


У(ху) -У(0.2) 2(1—0.4) — 0.832 
SWU =25kg (1.758 - 6.188) — 113.384 kg (0.832 - 6.188) 
= 496.535 kg 
(c) xp = 0.2, хт = 0.002, xp = 0.00711, Mp = 113.384 kg 


M, = 0 H3.384 kg = 4393.35 kg 
V(xg) = V(0.00711) = (1 0.01422) —=4869 


SWU = 113.384 kg (0.832 - 6.188) — 4393.35 kg (4.869 - 6.188) 
= 5187.54 kg 


31. (a) Derive an explicit expression for the mass of tails produced in a specified enrichment requiring a given 
amount of separative work. (b) Show that in enriching natural uranium to 3 w/o, approximately 1 kg of 0.2 w/o 
tails are produced per SWU expended. (c) In 1981, the total enrichment capacity in operation in the world was 
28.4 10° SWU/yr. At what rate was depleted uranium being produced, assuming all the capacity was used to 
produce 3 w/o fuel? 


[sol] (a) Mr = ОМт- 
Mp = Mp - Mr 
SWU = МьУ(хь) + МтУ(хт) - MFV (xr) <Eq. 4.36» 
= (My - Mr) V(xp) + МтУ(хт) - MFV (Xr) 
= М[У(хт) - V(xp)] - Mr[V Xr) — V(xp)] 
= М У(хт) – V(xp)] - OMHEV(xr) — У(хр)] 


Мт = _ 


(b) xy = 0.00711, хь = 0.03, хт = 0.002 
V(0.03) = 3.268 
V(0.002) = 6.188 
V(0.00711) = 4.869 


=1.04 = 1 
Approximately 1 kg of 0.2% tails are produced per SWU 
(с) Depleted uranium production rate is 28.4 107 kg/yr, because approximately 1 kg of 0.2% tails 
are produced per SWU. 


32. In 1980, the United States had approximately 300,000 tons of depleted uranium in storage at its gaseous 
diffusion plants. If the entire 1980 electrical capacity of 600 GWe were furnished by LMFBRs fueled with this 
uranium, how long would the 1980 depleted uranium stock (which is continuing to grow; see preceding problem) 
last? Use LMFBR data from Fig. 4.40. 


[sol] МЕ = 2,513 | —Mass flows in kg’s per 0.75 GWe-yr> 
= 2,010,400 kg 
xy = 0.00711, xp = 0.03, хт = 0.002 
Mp = OMr= ( — ) —2,019,400 kg = 366,898 kg 
Мт = My - Mp = 2,010,400 kg - 366,898 kg = 1,643,502 kg 
Lasting time = —=+82-537years 


33. Show that approximately 116,000 SWU of separative work is required annually to maintain the nominal 
1,000 MWe LWR described in Fig. 4.37. 


[sol] SWU = Mp[V(xp) — У(хт)] - МЕ[У(хк) — V(x1)] 
Mp = 27,249 kg, My = 149,297 kg 
хь = 0.03, хт = 0.002, хь = 0.00711 
V(0.03) = (1—0.06) | —3.268 


V(0.002) 2 6.188 
V(0.00711) = 4.869 
SWU = 27,249 Ко (3.268 - 6.188) – 149,297 kg (4.869 - 6.188) 


= 117356 kg = 116,000 kg 


34. The enrichment consortium EURODIF operates gaseous diffusion plants in France with a capacity of 10.8 
million SWU/year. Using the result of the preceding problem, how much LWR capacity can this company 
service? 


[sol] 1,000 MWe : 116,000 SWU/yr = x : 10,800,000 SWU/yr 
x = ———————— — —2-931093.4-MWe 


35. Compare and plot the relative cost in SWU of enriched uranium, per unit of contained ^ `U, as a function of 
enrichment to 93 w/o. Compare with Fig. 4.45. 


[sol] хт = 0.002, хр = 0.00711 
SWU = Mp[V(xp) — V(xy)] - Me[V(xr) — V(x1)] 
—z[V(xp) - Vxy]- (хв) - У(хү)] 
= [V(xp) — V(0.002)] - OEV«9.00711) — V(0.002)] 


=(1-2xp)  —— 6.188- 44.869— 6:188) 
-(1-2xp) —+  — 6488 

If xp = 0.93 

(xp = 0.93) = (1-2 0.93) —+ ——-6488 


= 235.83 
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36. Based on the projections of future installed nuclear electric generating capacity given in Table 4.8, how large 
an industry (dollars per year) can isotope enrichment be expected to be in 2010 based on current U.S. SWU 
prices if most of the capacity is in LWRs? 


[sol] 


37. The principle of enrichment by gaseous diffusion can be seen in the following way. Consider a chamber 
(diffuser) split into two volumes, A and B, by a porous barrier with openings having a total area 5 as shown in 
Fig. 4.48. Volume A contains two isotopic species of mass and at atom densities and , 

respectively. Volume B is at low pressure and is pumped out at the rate of F m3/sec. According to kinetic theory, 
the number of atoms or molecules in a gas striking the wall of a container per /sec 15 


where N is the atom density and is the average speed given by 





? 


where M is the mass of the atom and T 15 the absolute temperature. (a) Show that in equilibrium the atom 
concentration in B 18 


where апа refer to the two isotopic species. (b) Show that the ratio of equilibrium concentrations in B 15 
given by 


Number pumped out = — 


Number of atoms stibing the wall per second JS =(+)5 


38. The quantity 


in the preceding problem 18 called the ideal separation factor .Compute for a gas consisting of molecules 
of апа . (The small value of in this case means that a great many stages must be used to 
enrich uranium to useful levels.) 


[sol] 


39. If, using LIS, it is possible to enrich diffusion plant tails from 0.2 w/o to 3 w/o with residual tails of 0.08w/o, 
verify that the utilization of natural uranium for 3 w/o fuel would be increased by 18%. 


[sol] 
i) = == — ——0:1825- 
Mass of in = 0.00711 


Мавѕ оѓ іп = 0.030.1325=5.475 


ii) 
= —Ө=0-4825)—<0.033596 
= 0.2161 


Mass of іп =0.03 0.2161=6.4829 
Utilization = ——=914% 


40. To how much energy in e V does the isotopic shift of 0.1 shown in Fig. 4.52 correspond? 


[sol] 


2.46632731-2.466483669—0.000049062 


41. (a) Show that the specific activity of fuel irradiated for T days up to a burnup of B MWdlkg, t days after 
removal from a reactor, is given by 


[sol] 

Fission product activity = 1.4 [ — Ci.....(3.49) 
Burnup rate = 1.05P g/day...(3.57) Р g/day 
The fission of 1g yields 1 Муа 


В = — = 


PMw-B  MWd/kg- 


Mass of t operate on PMW for T days 


=(P = (Р 


(b) Compute and plot ће specific activity for fuel irradiated for 3 years to a burnup of 33MWd/kg from 1 day to 
] year after removal from the reactor 


[sol] 
B= 33MWd/kg, T=. 


i) t =Iday 


ii) t=365day 
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42. The fission yields of and are 0.0593 atoms per fission, respectively. Calculate the specific 


activity due to these nuclides and to in spent fuel irradiated to 33MWd/kg from the time the fuel is 
removed from the reactor up to 1,000 years. 


[sol] 
Fission rate 22.7 


For 
N=2.7 330.0593-5.2836 /kg 


1.5 Ci/kg 


For 
N=2.7 330.0623=5.551 /kg 


= 1.5 /kg 


For 


43. Carbon-14, with a half-life of 5,730 years, is produced in LWRs by way of an (n,p) reaction with nitrogen 
impurities (in both the fuel and the coolant water) and via an (n,) reaction with — . About 60 to 70 Ci of 

are generated this way per GWeyear. (a) What is the total activity in the fuel unloaded each year from a 

1 ,000MWe reactor that has operated with а 0.70 capacity factor? Assume that one third of the core is unloaded 
each year. (b) How much will be produced per year in the year 2010 by all the reactors in the world if high 
or low projections of Table 4.8 come to pass and all reactors are LWRs? 


(a) 

[sol] 0.7 capacity factor and 1/3 of the core is unloaded, the total activity is 
0.7 = 

b) 


[sol] 

assuming that about 65Ci of опе generated рег GWe 
i) for low projections 

activity — 65 


N= = 
M==12390g=12.39kg 





ii) for high projections 


activity = 65 


N- 


M==17.49k¢ 





Chapter 5 


Neutron Diffusion 
and Moderation 


1. A point source emits S neutron/sec isotropically in an finite vacuum. 
(a) Show that the neutron flux at the distance r from the source is given by 


[sol] It is important to note that this solution is quite different than would be obtained for a point source 
in a vacuum. This difference occurs since the diffusing media both absorbs and scatters neutrons, whereas 
a vacuum does not. 


So, — 


(b) What is the neutron current density vector at the same point? [Note: Neutrons do not diffuse in a vacuum. | 
[sol] - —— -- 


2. Three isotopic neutron source, each emitting S neutons/sec, are located in an infinite vacuum at the three 
corners of an equilateral triangle of side a. Find the flux and current at the midpoint of one side. 


[sol] reference problem (Example 5.3) 
The flux 2/а 


The current, 


(Due to characteristics of vector, other two currents are canceling each other.) 


3. Using Eqs. (5.10) and (5.11), estimate the diffusion coefficients of (a) beryllium, (b) graphite, for 
monoenergetic 0.0253 eV neutrons. 


[sol] Equ 5. 10 and 11 : 


(a) Beryllium 


(b) Graphite 


4. The neutron flux in a bare spherical reactor of spherical reactor of radius 50 cm is given by 


where r is measured from the center of the reactor. The diffusion coefficient for the system is 0.80 cm. (a) What 
is the maximum value of the flux in the reactor? (b) Calculate the neutron current density as a function of 
position in the reactor. (c) How many neutrons escape from the reactor per second? 


[sol] (a) 


(b) 


(c) Leakage rate 


5. Isotope point sources each emitting S neutrons/sec are placed in an infinite moderator at the four corners of a 
square of side a. Compute the flux and current at the midpoint of any side of the square and at its center. 


a 
[sol] —— S1 52 
= - P2 
(a) 2 —— --- = — О P1 
= E -- S3 
i JA J1 
шигээ [ м. 
P1 
(b) at P1 
at P2 


6. Find expressions for the flux and current at the point P2 in Fig. 5.5 


[sol] 
(a) ——— —=— 


(b) 


7. An isotropic point source emits S neutrons/sec in a n infinite moderator. 
(a) Compute the net number of neutrons passing per second through a spherical surface of radius r centered 
on the source. 
[sol] 
—.. (Eq 5.33) 


net number of neutrons passing per second through a spherical surface of radius r centered 


(b) Compute the number of neutrons absorbed per second within the sphere. 


[sol] 


(c) Verify the equation of continuity for the volume within the sphere. 


[sol] 


8. Two infinite planar sources each emitting S neutrons/cm2 are placed parallel to one another in an infinite 
moderator at the distance a apart. Calculate the flux and current as a function of distance from a plane midway 
between the two. 


[sol] 
i) flux is scalar factor. So the flux at x =0 


ii) from the planar source, current vetors have same magnitude and opposite direction. So current at x = 0 
is zero. 


Хо) = 0 


9. Suppose the two planar sources in the preceding problem are placed at right angles to one another. Derive 
expressions for the flux and current as a function of distance from the line of intersection of the sources in a 


plane bisecting the angle between the sources. 


[sol] 
Flux from the distance r 


? 


10 An infinite moderator contains unifonnly distributed isotropic sources emitting S neutrons/cm3 -sec. 
Detennine the steady-state flux and current at any point in the medium. 


[sol] 


Netron source in dv = Sdv 


dv 


‚ J(P)-0 


11. An infinite bare slab of moderator of thickness 2a contains unifonnly distributed sources emitting S 
neutrons/cm3 -sec. (a) Show that the flux in the slab is given by 





where x is measured from the center of the slab. (b) Verify the equation of continuity by computing per unit 
area of the slab the total number of neutrons (1) produced per sec within the slab; (11) absorbed per second within 
the slab; and (11) escaping per second from the slab. [Hint: The solution to an inhomogeneous differential 
equation is the sum of solutions to the homogeneous equation plus a particular solution. Try a constant for the 
particular solution.] 


[sol] 
(a) 


At slab geometry 


i) homogeneous solution 


ii) particular solution 


0- — = = 


boundary condition 


i) 


Sinhx = -sinh(-x), conshx = cosh(-x), A=0 


il) 


(b) 


i) S neutron / -sec 


ii) 


iii) J(x) = -D— —- — 





DivJ = D— S 


Continuity Eq 


12. A point source emitting S neutrons/sec is placed at the center of a sphere of 
moderator of radius R. (a) Show that the flux in the sphere is given by 


where r is the distance from the source. 


[sol] 
No source with r 


? 


Let w=r 


= C 
C B 
2 _ 2 шог р 
В.С 
EROS E 
- €D) = -4 - 

EN ES F 
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(b) Show that the number of neutrons leaking per second from the surface of the sphere is given by 


No.Leaking/sec = 


[sol] 
J= -D— | — — —— Ov ----- - - 


d/L 
cosh(d/L)  ,sinh(d/L) 


(c) What 1s the probability that a neutron emitted by the source escapes from the surface? 


[sol] 


13. An infinite planar source emitting S neutrons/cm^-sec is placed between infinite slabs of beryllium and 
graphite of thickness a and b, respectively, as shown in Fig. 5.9. Derive an expression for the neutron flux in the 
system. [Note: Since the media are different on opposite sides of the source, this problem is not symmetric and 
the source condition Eq. (5.28) is not valid. The appropriate boundary conditions for this problem are 


(i) lim, [Øx > 0) -Øx < 0)] 20 
(ii) lim, ох > 0) - x < 0)] =S 


Condition (1) states, in effect, that 0 is continuous at the source, whereas (ii) accounts for the neutrons emitted 
from the source.] 


Planar source emitting S 
/ neutrons/cm*- sec at x = 0 
x 


Air Beryllium Graphite Air 


Figure 5.9 Infinite planar source 
between two finite slabs of material. 


[sol] — Ө=0,х #0 «Eq. 5.25» 
general solutions : O4 < = C4cosh—- C,sinh — 
OQ. о = C3cosh—+ Cysinh — 
Let, =a +d, =b+d, 
Boundary Condition 1: 0,0 = 0, Ø )=0 
ØQ = C,cosh — C;sinh —=0 


C; = С; 





@,() = Cacosh — Cysinh — 0 


С. = -C,—— 


Oo = Ci(cosh- —sinh }- 


Өк» о = C3(cosh— — nh 3} 
Boundary Condition 2 : lim, ,4[O0(x > 0) – О(х < 0)] =0 


lim, .,o[Cs( cosh— — sinh j— C,(cosh —+ — sinh )]20 
C3(1 - 0) —C,(1 — 0) = 0, Сз = Cj 


Øx <0) = Ci(cosh— — sinh у- 
О.о = Ci(cosh—- —sinh >} 


Boundary Condition 3 : і, ,([Ј(х > 0) – Хх < 0)] =S 


lim,.o[-D:Csinh — ——eosh }+D,C, (sinh — —eosh }}=S 


D;,C,— ——- D;C; pcc S 

C, = 

Solution : 

Ox <0) = —————— — — (cosh ++ — sinh r 


шэг 115 


9, >0) = (cosh — —sinh Э- 


14. A sphere of moderator of radius R contains uniformly distributed sources emitting S neutrons/cm -sec. (а) 
show that the flux in the sphere is given by 


(b) Derive an expression for the current density at any point in the sphere. (c) How many neutrons leak from the 
sphere per second? (d) What is the average probability that a source neutron will escape from the sphere? 


[sol] (a) —éà — Ø=-  -«Eq. 5.19, 5.30» 


Let, w =r@ 
—— 


(1) if homogeneous, 
W, = Asinh-+ Bcosh - 


QO, = -sinh+ osh - 


(ii) if particular 


W, = Cr 

— ww. = 
(ewe = 
Ww, = — 
Ø, = — 

Ø = Øn +0, 


= -sinh+ eosh+  — 
Boundary Condition 1: = finite 


0-0--4 —— finite 

В=0 

Ø= sinh+ -- 

Boundary Condition 2: Oa) = 0 

smh --4 —- 

ав. 

0--8------42- > 
= Шш m — 


(b) current density J «Eq. 5.8» 
--DO--D t 4- — ——3M4 


--Di——Cr-2 — ——r-1. ———cosb} - 





—— — 


(c) leakage rate 
= 4R?J 


= 4R? ——Ó ——.  —Jlumeutrons/s 


(d) The probability that a neutron escapes is equal to the number that escape per cm/sec divided 
by the number emitted рег ст2/ѕес by the source 





Е > 


The number that escape рег cm?/sec = 4R? 


The number emitted per cm?/sec by the source = S 
The probability that a neutron escapes = 


= 4R? —— — —- 


15. The three-group fluxes for a base spherical fast reactor of radius К = 50 cm are given by the following 
expressions: 








Qr) = sin ( — 
Or) = sin ( — 
Qr) = sin (— 





The group-diffusion coefficients are D; = 2.2 cm, D» = 1.7 cm, and D; = 1.05 cm. Calculate the total leakage of 


neutrons from the reactor in all three groups. [Note: Ignore the extrapolation distance.] 


[sol] leakage in group е: div:J 
Total leakage of neutrons : = - D,’O,— D//0,— D7 Ø; 
20, = —t sin ( + 
-(3 10^) ба’ >} 
=(3 10°) Gr’sin —+г! eos Э 
-3 10^ Отт  — r° eos — r° eos —4r! -Csin 9) 





«Ignore the extrapolation distance? 
-3 10° ( + y — — 


“0, = -0,-0--- 





1.05 — 








Total leakage of neutrons = -2.2 1.7 
--2.57 10" neutrons 


16. The thermal flux in the center of a beam tube of a certain reactor is 2 10 neutrons/cm?-sec. The 
temperature in this region is 150°C. Calculate (a) the thermal neutrons density; (b) the energy Ет; (c) the 2,200 


meters-per-second flux. 


[sol] (a) От = nx, «Eq. 5.56» 


n= — 
vr 21.284T'7 10“ cm/sec «Eq. 5.55» 

= 1.284(273 + 150)'? 10“ cm/sec 

=2.64 10° cm/sec 
n = (2 х 107 neutrons/em-sec x )7(2 2.64 10° cm/sec) 


=6.71 10’ neutrons/cm? 


(b) Ey =8.617T 10° eV «Eq. 5.54» 
= 8.617(273 + 150) 10° eV 


= 0.0345 eV 


(с) О, = nv, <Page 253» 
= 6.71 10’ neutrons/em*? 2.2 10° cm/sec 
=1.48 10" neutrons/cm?-sec 


17. The thermal flux at the center of a graphite research reactor is 5 10” neutrons/cm’-sec. The temperature 
of the system at this point is 120°C. Compare the neutron density at this point with the atom density of the 


graphite. 


[sol] От = WT «Eq. 5.56» 


n= — 
ут = 1.2841"? 10“ cm/sec «Eq. 5.55» 

= 1.284(273 + 120)'* 10“ cm/sec 

=2.55 10 cm/sec 
nz(5x 1077 neutrons/cm^-sec x )/ (2 2.55 10° cm/sec) 


= 1.74 10’ neutrons/cm? 
atom density of the graphite = (1.6 g/cm3 0.6022 107 atoms) / 12g «1.6 g/cm? – Page 254» 


= 8.03 10” atoms/cm? 


18. The thermal flux in a bare cubical reactor 18 given approximately by the function 


@т(х, у, 2) = A cos C cos ( cos ( 7 


where A is a constant, a is the length of a side of the cube, à is a plus 2d, d is the extrapolation distance, and x, y, 
and z are measured from the center of the reactor. Derive expressions for the (a) thermal neutron current as a 
function of position in the reactor; (b) number of thermal neutrons leaking per second from each side of the 
reactor; and (c) total number of thermal neutrons leaking per second from the reactor. 


[sol] (a) thermal neutron current J 
--00) 
3 (х,у,?) = -DG+ i t )O(x,y,z) 

= -D@(x,y,z)i - D -O(x,y,z)j - D O(x,y,z)k 

= -DiAcos( —cos( —cos( Н 
- DfAcos( cos( Эсоѕ( Jj 
- DfAcos( —cos( —cos( ]k 

= DAdsin( cos( —cos( і + cos( sin( cos( -} + cos( cos( sin( ЭК) 


(b) each six surface is same. 
number of thermal neutrons leaking per second from each side of the reactor 


= == — _ _ —dz 
Bob = = = - 
-DA x 2 х-х 2= —neutrons/sec 


(c) total number of thermal neutrons leaking per second from the reactor 
neutrons/sec 





= 6 Х өөө: 
19. A planar source at the center of an infinite slab of graphite 2 meters thick emits 10? thermal neutrons per 
cm’/sec. Given that the system is at room temperature, calculate the: (a) total number of thermal neutrons in the 
slab per ст? at any time; (b) number of thermal neutrons absorbed per cm/sec of the slab; (c) neutron current as 


a function of position in the slab; (d) total number of neutrons leaking per cm /sec from the two surfaces of the 
slab; (e) probability that a source neutron does not leak from the slab. 


[sol] (a) 


-2AL = — 


(d) А 


(e) Probability —— 


20. The thermal flux in a bare cubical reactor of side — is given by the expression 


ü 2 = ЭР 
where x,y, and z are measured from the center of the reactor апа . The temperature is 400°C, and the 
measured values of the thermal diffusion coefficient and diffusion length are 0.84cm and 17.5cm, respectively. 
(a) How many moles of thermal neutrons are there in the entire reactor? (b) Calculate the neutron current 
density vector as a function of position in the reactor. (c) How many thermal neutrons leak from the reactor per 
second? (d) How many thermal neutrons are absorbed in the reactor per second? (e) What is the relative 
probability that a thermal neutron will leak from the reactor? 


[sol] (a) 


(Total 8) — — — 


(b) = = = == = 


(c) Leakage rate 


(d) Absorption rate — - I 


(e) — —— 


21. The thermal (0.0253 eV) cross-section for the production of a 54-minute isomer of ''°In via the y 
reaction with "Эш is 157b. A thin indium foil weighing 0.15g is placed in the beam tube described in Problem 
5.16. (a) At what rate are thermal neutrons absorbed by the foil? (b) 

After 1 hour in the beam tube, what is the 54-minute activity of "бп? [Note: 115In has an isotopic abundance of 
95.7% and is a non-1/v absorber. | 


[sol] (Table II.3) 
(Table 3.2) ° 


(a) = ын 


(9) 


22. The thermal flux in a bare spherical reactor 1m in diameter is given approximately by 
— ———4geutrons/cm-sec. 


If the reactor is moderated and cooled by unit density water that takes up to one-third of the reactor volume, how 
many grams of ^H are produced per year in the reactor? Assume that the water is only slightly warmed by the 
heat of the reactor. 


[sol] Table II.3 


(Absorption rate) = = 


(Absorption rate) 


23. Tritium (ЗН) is produced in nuclear reactors by the absorption of thermal neutrons by 14 via the reaction 
“Г(а,а)7 H. The cross-section for this reaction at 0.0253eV is 940b and it is l/v. (a) Show that the annual 
production of 3H in a thermal flux per gram of °Li, is given by 


Production rate/g= — atoms/g-yr. 

(b) Compute the annual production in curies of ^H per gram of 614 in a flux of neutrons/cm"- 
sec. 

[sol] 


(a) Production rate-absorption rate of thermal neutrons by 


= —=1.22 


(b) 


Production- 


24. A radioactive sample with half-life 1$ placed in a thermal reactor at a point where the thermal flux is 


. Show that the sample disappears as the result of its own decay and by neutron absorption with an effective 
half-life given by 
б 


where o isthe average thermal absorption cross-section of the sample. 


[sol] 


25. Calculate the thermal diffusion coefficient and diffusion length of water near the outlet of a pressurized- 


water reactor, that is at about 300°C and density of 0.68g/cm”. 


[sol] 
(р, T,)=0.l6em, p,=1g/cm’, p-0.68g/cm? 
T= 300+273=573°K, T,=293°K, т-=0.470 
L^; (Po, T,)=8. lem” 
i) (р, = (Po, T)(pJp)(T/ Т.)" 

= 0.16 x (1/0.68)(533/293)°479=0.312cm 
ii) L^7 (р, T)= L^; (Po, То) (Po/P)2(T/ То)" +12 


=8.1 x (1/0.68) 7 (533/293)??? 
-23.2cm? 


Lr (р, T)24.8cm 


26. Calculate the for natural uranium at room temperature and 350°C the value of >,, D, and Ly. The measured 
value of Іт at room temperature is 1.55 cm. [Note: The value of os is constant at low neutron energies and 18 


approximately the same for 235U and 238U. The densities of U at 350°С and room temperature are essentially 


the same] 


[sol] 
P= р,, T,2293K ,Т-623К, m=0, 6,(235U)=680.8b, 6,(238U)=2.70b 
У.-0.4301сш" 


2.(Т)-+ i) for 2350, 0.9780 at T, , 0.9335 at T 


ii) for 239U, 1.0017 at T, , 1.0067 at T 
6,(235U)=680.8b .,0,(238U)-2.70b 


i) (Po, T2 4/3 =1/3У,(1- 
-1/3 x 0.4301 x (1-2.801 x 107), where -=2/3 x 238-1-2.801 x 107 
-20.777cm'! 
ii) L^; (р, Т,)=1.552=2.4025ст? 
L^; (p, Т)= )= L^; (Po, То) (p/p)^(T/ T,)"*!? 
= 2.4025 х (693/293)1/2=3.5ст* , 1л-1.87сш 
iii) (ТЕХ [g,(T) У, (T,/ T) 7] | N(U)20.04833 x 10”#/cm3 
= (8 ”,М(2350) 6,(235U)+ g% М№(2380) с,(2380) 
- (0.9780 x 0.0072 х 0.04833 х 680.8-1.0017 х 0.9927 х 0.04833 х 2.70) 
= 0.32 cm” 
(Т)= ga(To) У„(Ео) (Т,/Т) У? 
-[(g7?,N(235U) в. (2350)+ 27 М(2380) 6,(238U)] x (293/623) 12 
=[(0.9780 x 0.0072 x 0.04833 x 680.8+1.0017 х 0.9927 x 0.04833 x 2.70)] x (293/623) 1? 
20.215 cm” 
(Т,)= (Pos T)/ L^; (Pos T,)20.777/1.55?20.32cm'! 


(T)= (p, T)/ L^, (p, T)20.777/3.520.22cm ! 
27. Repeat the calculation of problem 5.26 for U enriched to 2 weight/percent in 235U. 


[sol] 
Volume percent of 235U 


= [(0.02 x 235.0439-0.98 x 238.0508)/235.0439] x 0.02=0.0203 
=2.03percent 


6,(235U)= 6,(238U)= o, at low energy 
$.- М(2350) 6,(235U)+ N(238U) 6,(238U) 
= [N(235U) +N(238U)] х бс, 
=(19.1 х 0.602217 x 8.9)/0.01(2.03 x 235.0439497.97 x 238.0508)=0.43015 cm" 


i) (Pos То)= №/3 21/35 (1) 
=1/3 x 0.43015 x (1-2.801 х 107), where -—2/3 x 238=1-2.801 x 10° 
=0.777 cm” 


ii) L^; (р, Т,)=1.552=2.4025ст? 
L^t (p, T) = Ly (Pos То) (Po/p)2(T/ Т)" 
= 2.4025 х (693/293)1/2=3.5ст* , 1л-1.87сш 
iii)  (T)-g,(T) Ya(Eo) (T, / T) V? N(U)=0.04833 х 10°°#/ст3 


= (8 ”,М(2350) 6,(235U)+ g% М№(2380) с,(2380) 
- (0.9780 x 0.0072 х 0.04833 х 680.8-1.0017 х 0.9927 х 0.04833 х 2.70) 
= 0.32 cm” 
(Т)= ga(To) Y,(E) (Т,/Т) !? 
-[(g7?,N(235U) 6,(235U)+ 2? М(2380) 6,(238U)] x (293/623) 12 
=[(0.9780 x 0.0072 x 0.04833 x 680.8+1.0017 x 0.9927 x 0.04833 x 2.70)] x (293/623) 1/2 
-0.215 cm” 
(Т,)= (Po, T)/ L^; (Pos T,)20.777/1.55?20.32cm'! 
(Т)= (p, T)/ L^, (p, T)20.777/3.5-0.22cm ! 


28. Compute and plot for mixtures of D;O and H,O, at H5O concentrations up to 5 weight/percent: (a) У, at 
0.0523 eV; (b) L4 at room temoerature. 


[sol] 
Let weight percent of H,O x, volume percent of H;O a 


M=1/100[18.0153a+20.0276(100-a)]={[(x/18.0153)+(100-x)/(20.0276)] x 1/100]! 
N=pN,/M=(20.0276 x 0.6022 x 107/[107(2002.76-2.01232a)]2(12.06 x 107*)/(20.0276-0.020123a) 


a={[(x/18.0153)+(100-x)/(20.0276)] x 1/100] x (1/18.0153)x 


= 995.259-(8.91 x 10°/x+895.26) 
i) Y,2 N(D;O) 6a(D20)+ N(H2O) б,(Н:О) 
=1/100[(100-a)N o,(D,O)+ a)+aN б,(Н:0)| 


=1/(2002.76-2.0123a) х [(100-a)0.00133+0.644a] 
=(0.329a+0.066)/(a-995.259) а-995.259-(8.91 х 10°/x+895.26) 
-(327.442--5.189)/8.91 х 10” 
-3.675 х 10 х--5.824 x 10° 
1) (D;O) -0.87сш 
L'r-/Ya- /хв (1) XE) (T/T) 17 
=0.87/[ х 1 х (3.675 х 10 х--5.824 x 109) x1 
-0.87/(3.257 х 10^x--4.683 х 10°) 
-2.671 х 10°/(x+0.0144) 


Ly =[2.671 х 10/(x40.0144)]'? 


29. Calculate the thermal neutron diffusion length at room temperature in water solutions of boric 
acid(H3BO3)at the following concentrations: (a) 1Og/liter, (b) Ig/liter and (c) 0.1 g/liter. [Hint: Because of the 
small concentration of the boric acid, the diffusion coefficient for the mixture is essentially the same as that of 
pure water. | 


[sol] 
(НО) =1g/cm?, =0.16ст 
Let weight percent of НВО; x, volume percent of HBO; a 


M=1/100[18.0153(100-a)+(3-1.00797+10.811+3-15.9994)a] 
= 0.438а+18.0153 
= [(100-x)/(18.0153)+ (x/(3-1.00797+10.811+3-15.9994)]" 


N=pN,/M=(0.602210™)/ (0.438a+18.0153) molecules/cm? 


a=[(100-x)/(18.0153)+ (x/(61.833)] -(x/61.833) 
= 100x/(343.2-2.433x) 


i) 10g/ cm? = 0.01g/cm? — 1 w/o 
а-100/(343.2-2.433)-0.293 
с.(В) o,( НВО) = 759b 
> а= N(H50) в.(Н2О)+ М(НзВО:) ба (НэВОЗ) 
=[(100-0.293)N б,(Н:О) + 0.293 М c,(H3BO3)]-1/100 
Х-4 х 0.6022 x 10?)/ (0.438 x 0.293+18.0153)=0.0332 х 10” 


У,410100-0.293):0.0332-0.664 + 0.293-0.0332- 7591-1/100 
-20.0958cm'! 


L^, =/У»=/[ X г„(Т) Ya(Eo) (ТТ) ?] 
=0.166/(-0.0958) 


-1.884cm? 
Lr = 1.372cm 


ii) 1g/ cm? = 0.001g/cm? — 0.1 w/o 
a=100-0.1/(343.2-2.433-0.1)=0.029 w/o 
N=(0.6022 х 10”4)/ (0.438 x 0.029+18.0153)=0.0334 х 1074 


У,4111000-0.029)-0.0334-0.664 + 0.029-0.0334- 759]-1/100 
-0.0295cm'! 


L^, =0.166/(-0.0295) 
-26.12cm? 
Гл = 2.47ст 


iii) 0.1g/ cm? = 0.0001g/cm^ — 0.01 w/o 
a=100-0.01/(343.2-2.433-0.01)=0.0029 w/o 
N=(0.6022 х 1029/ (0.438 х 0.0029+18.0153)=0.0334 x 10” 


У,41100-0.0029)-0.0334-0.664 + 0.0029-0.0334- 759]-1/100 
=0.0229em! 


L^, =0.166/(-0.0229) 
=7.884ст? 
Lr = 2.8cm 


30. An infinite slab of ordinary water 16 cm thick contains a planar source at its center emitting 10* thermal 
neutrons рег cm’/sec. Compute and plot the thermal flux within the slab. 


[sol] —————————— — 
Where) S=10°cm7/sec, Їлг-2.85сш, =0.16cm, a-8cm, 4-2.13 =0.341ст 
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31. Repeat the calculations of Problem 5.30 and plot the thermal flux for boric acid solutions having the 
concentrations given in Problem 5.29. 


[sol] at L;=1.884cm 


at L7=2.47cm 


at L;=2.8cm 
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32. Verify that Eq. (5.72) is in fact the solution to Eq. (5.68). Discuss this solution in the cases where тг > 0 
and tr > LA 


[sol] Eq.(5.72), —)- 
Eq.(5.68,—— — — 


attr — 0, _— 
33. Calculate the neutron age of fission neutrons in water for the conditions given in Problem 5.25. 


[sol] — —m 


34. An infinite slab of moderator of extrapolated thickness à contains sources of fast neutrons distributed 
according to the function 


Using the two-group method, derive an expression for the thermal flux in the slab. 


[sol] — 


= = --(1) 


Particular Solution 


= — .-->insert the above eq(1) 


- = = = = E — 
i= 
{ - = 


———eos( — 


Homogeneous solution 


— — ———eoes( —=А“ ---608( — 
Boundary condition 
1) - - 
A 
А“ i i —— —eos( — 
2) 
= = и ш ш - ——eos( — 
=> = i i 


Chapter 6 


Neutron Reactor Theory 





1. Calculate the fuel utilization and infinite multiplication factor for a fast reactor consisting of a mixture of 
liquid sodium and plutonium, in which the plutonium is present to 3.0 w/o. The density of the mixture is 


approximately 1.0 g/cm3 


[sol] — — 


2. The core of a certain fast reactor consists of an array of uranium fuel elements immersed in liquid sodium. 
The uranium is enriched to 25.6 w/o in 235 U and comprises 3796 of the core volume. Calculate for this core (a) 


the average atom densities of sodium, 25U. апа ?5U; (b) the fuel utilization; (c) the value of n; (d) the infinite 


multiplication factor. 


[sol] (a) 


(d) = 


3. A bare-cylinder reactor of height 100 cm and diameter 100 cm is operating at a steady-state power of 20 MW. 
If the origin is taken at the center of the reactor, what is the power density at the point r = 7 cm, Z = -22.7 cm? 


[sol] — —— = 


4. In a spherical reactor of radius 45 cm, the fission rate density is measured as 2.5 x 10" fissions/cm"-sec at a 


point 35 cm from the center of the reactor. (a) At what steady-state power is the reactor operating? (b) What is 
the fission rate density at the center of the reactor? 


[sol] (a) 


(b) р 


5. Using the flux function given in Table 6.2 for a critical finite cylindrical reactor, derive the value of ће 
constant A. 


[sol] — 2 


6. The core of a certain reflected reactor consists of a cylinder 10-ft high and 10 ft in diameter. The measured 
maximum-to-average flux is 1.5. When the reactor is operated at a power level of 825 MW, what is the 
maximum power density in the reactor in kW/liter? 


[sol] === 


7. Suppose the reactor described in Example 6.3 was operated at a thermal power level of 1 kilowatt. How many 
neutrons would escape from the reactor per second? [Hint: See Example 6.4.] 


[sol] = 


Non-leakage probability 


Leakage probability 


Number of produced neutrons by fission reaction 





9. Show that in a one-group model, the power produced by a reactor per unit mass of fissile material is given by 


Where is the one-group fission cross-section is the average one-group flux, is Avogadro's number, 
and 15 the gram atomic weight of the fuel. 


[sol] A one-group model 


9. (a) Estimate the critical radius of a hypothetical bare spherical reactor having the same composition as the 
reactor in Problem 6.1. (b) If the reactor operates at a thermal power level of 500 MW, what is the maximum 
value of the flux? (c) What is the probability that a fission neutron will escape from the reactor? 


[sol] Table II.3 
Table 6.1 





D=2.13D=8.462cm 


(Probability)=1-0.432=0.568 


10. An infinite slab of moderator of thickness 2a contains at its center a thin sheet of U of thickness t. Show 
that in one-group theory the condition for criticality of this system can be written approximately as 


Where D and L are the diffusion parameters of the moderator and is te macroscopic absorption cross- 
secton of the 


[sol] «Fuel» 


«Moderator» 


The condition for criticality of this system 


11. A large research reactor consists of a cubical array of natural uranium rods in a graphite moderator. The 
reactor is 25ft on a side and operates at a power of 20MW. The average value of of (а) 

Calculate the buckling. (b) What is the maximum value of the thermal flux? (c) What is the average value of the 
thermal flux? (d) At what rate is 770 being consumed in the reactor? 


[sol] (a) - Е 7 — — — 
(b) 


0) - : = " = 
(4) 


12. Show with a recoverable energy per fission of 200MeV the power of a 235U-fueled reactor operating at the 
temperature T given by either of the following expressions : 


Or 
where is the total amount in kg of 235U in the reactor, is the non-1/ factor for fission, is the 
average 2,200 meters-per-second flux, and is the average thermal flux. 


[sol] Table 3.4 


13. Solve the following equations: 


[sol] Dz 


? 


14. Solve the following equations: 


[sol] D=0 
(Trivial solution) : (x,y,z)=(0,0,0) 
x-5.0y=9.0z 
7х+4.5у=-8.12 


15. A homogeneous solution of 2350 and ; contains 10 grams of 2350 per liter of solution. Compute (а) 


the atom density of 235U and the molecular density of ; (b) the thermal utilization; (c) the thermal diffusion 
area and length; (d) the infinite multiplication factor. 


[sol] 

(a) 

N(U-235) = — — — ——-2.56x atoms/c 
М(О) =  — —— ——— molecules / c 


(b) 


(c) 
= (1-0.938) x 8.1 = 0.5с 


(d) 


16. Compute the thermal diffusion length for homogeneous mixtures of 235U and the following moderators at 
the given fuel concentrations and temperatures. Graphite: М (25) / М (e) = 4.7 x 1 0-6; T = 200DC. Beryllium: N 
(25) / М (Be) 1.3x 10-5 ; T= 1 00DC. D20: М (25)/N(D20) 1 .4 x 1 0-6; T = 20DC. H20: М (25) / М (H20) = 
9.2 x 1 0-4; T= 20DC. 


[sol] 
(a) 
Z= = — —————47x x 


(d) 2 = 9.2 x 


17. Consider а critical bare slab reactor 200 cm thick consisting of a homogeneous mixture of 2350 апа 
graphite. The maximum thermal flux is 5 X 1 012 neutrons/cm2-sec. Using modified one-group theory, calculate: 
(a) the buckling of the reactor; (b) the critical atomic concentration of uranium; (c) the thermal diffusion area; (d) 
the value of  ; (e) the thermal flux and current throughout the slab; (f) the thermal power produced рег cm2 

of this slab. 


[sol] 
(a) 


-=Z =— 
2 
(c) 
= —=— -=293x cm 

(d) з= —=2.056 x —=1,296 
(e) — —#/ 

Ja)--D- — - 
(D A-,—R- | —— 

+ = 82.2 х х0.9759=3.97х с 
Р= MW/c =8.1 W/c 


18. The binding energy of the last neutron in 13C is 4.95 Me V. Estimate the recoverable energy per fission in a 
large graphite-moderated, 235U-fueled reactor from which there 15 little or no leakage of neutrons or y-rays. 


[sol] 


Recoverable energy = 200-4.95 = 195.05 MeV = 3.12x Joules 


19. Calculate the concentrations in grams per liter of ( 1 ) 235U, (2) 233U, and (3) 239pu required 
for criticality of infinite homogeneous mixtures of these fuels and the following moderators: (a) H20, (b) D20, 


(c) Be, (d) graphite. 


[sol] 


Z2 2— — = —— 


(235) = 2.065 (233) = 2.284 (239) = 2.035 
7235) = 0.739, Z(233) = 0.779, 7,239) = 0.966 


Ga(235) = 0.9780, ga(233) = 0.9983, ga(239) = 1.0723 


For 
A= 0.939 x —-0:331— 


(a) 
—=1.22x gm/ 


(b) 


Ж — = 0.024 


x ——= 0.63 


x ——0.15 


For 
A= 0.779 x —-0:314— 


——= 11.6 


X —=).023 
x — —z 0.59 
X — 2-0.14 


For 
A= 0.966 x —-0:213— 


——= 7.86 


X — ——0.015 
X —= 0.4 


X — —-0.0 


20. A bare-spherical reactor 50 cm in radius is composed of a homogeneous mixture of 235U and beryllium. The 
reactor operates at a power level of 50 thermal kilowatts. Using modified one-group theory, compute: (a) the 
critical mass of 235U; (b) the thermal flux throughout the reactor; (c) the leakage of neutrons from the reactor; 
(d) the rate of consumption of 235U. 


[sol] 
(а) = - 


(b) 
at 20 


(c) 
J(R) = -D— 


=4 
(d) consumption rate 
= 1.05(1+)P gm/day = 1.05 x (1+0.169)x50x 











21. The flux in a bare-finite cylindrical reactor of radius г and height H is given by ¢ = Alo ((12.405r ) cos (rriiz ) 
Find A if the reactor 18 operated at a power of P watts. 





[sol] 
—=30x— —20.019429 m/ 
— — ——=ыЫ9х 
Z = =+49x 
= (1-0.628) x 8.1 = 3.013c 
i) For Spherical tank 


ii) For cylindrical tank 


22. It is proposed to store H2O solutions of fully enriched uranyl sulfate (235002504) with a concentration of 
30 g of this chemical per liter. Is this a safe procedure when using a tank of unspecified size? 


[sol] —Ó à p _ 


For spherical tank 


For cylindrical tank 


23. Show that the flux in a bare-cubical reactor of side a is 


[sol] 017|0| ==} Ало 


24. If the reactor in Problem 6.23 is operating at P watts, show that the constant A is 


[sol] 0j7|0]|-- 4] E ei eor SAI 


25. A bare-thermal reactor in the shape of a cube consists of a homogeneous mixture of 235U and graphite. The 
ratio of atom densities is Nf/Nm = 1.0*10-5 and the fuel temperature is 250°C. Using modified one-group 


theory, calculate: (a) the critical dimensions; (b) the critical mass; (c) the maximum thermal flux when the 
reactor operates at a power of 1 kW. 


[sol] 


26. The origin version of the Brookhaven Research reactor consisted of a cube of graphite that contained a 
regular array of natural uranium rods, each of which was located in an air channel through the graphite. When 
the reactor was operated at a thermal power level of 22 MW, the average fuel temperature was approximately 
300°C and the maximum thermal flux was 5*1012 neutrons/cm2-sec. The average values of апа were 

325 cm2 and 396 cm2, respectively, and = 1.0735. (a) calculate the critical dimensions of the reactor; (b) 
What was the total amount of natural uranium in the reactor? 


[sol] <= ЕЛУЗЧҸҮ 


27. Using one-group theory, derive expressions for the flux and the condition for criticality for the following 
reactors: (a) an infinite slab of thickness a, infinite reflector on both sides; (b) an infinite slab of thickness a, 
reflectors of thickness n on sides; (c) a sphere of radius R, reflector of thickness b. 


[sol] 


b) 
1) core region 


2) reflector region 


2) Interface condition 


с) 


1) core region 


2) reflector region 





3) interface condition 


28. The core of a spherical reactor consists of a homogeneous mixture of 235U and graphite with a fuel- 
moderator atom ration Nf/Nm = 6.8 * 10-6 The core is surrounded by an infinite graphite reflector. The reactor 


operates at a thermal power of 100 kW. Calculate the: (a) value of koo; (b) critical core radius; (c) critical mass; 
(d) reflector savings; (e)thermal flux throughout the reactor; (f) maximum-to-average flux ratio. 


[sol] 


29. Estimate the new critical radius and critical mass of a reactor with the same composition as described in 
problem 6.20 when the core is surrounded by an infinite beryllium reflector. 


[Sol] 

рк=шь/У=1741/—660)5=3.325 x 10?g/cm? 
7=4.98 

f=Z/Z+1=0.083 

К =e7f =2.065 x 0.833=1.73 


Ly2=(1-f) Lyy’?=(1-0.833)480=80. 16cm” 
B?-k -1/ Lqzc270.72/80.16—8.892-10?-cm? 
B-20.0948cm'! 


Sjince the moderater is the same as reflector, 


BcotBR=-1/ Їл, 
cotBR--1/0.0948-21--0.502 
ВК=2.23 
R=2.23/B=2.23/0.0948=23.5cm 
Critical mass 
mp=—23.5cm)*-3.325-10°g/cm° 
-181.7g 


30. Show that the solution to the one-speed diffusion equation in an infinite cylinder made of a nonmultiplying 
media as 


Where Ko and J1 are modified Bessel functions. 


3]. An infinite-cylinder reactor core of radius К 15 surrounded by an infinitely thick reflector. Using one-group 
theory, (a) find the expressions for the fluxes in the core and reflector; (b) show that the condition for criticality 
is 


DcB x {J1(BR)/Jo(BR)= Dr/Lr x {K1(R/Lr)/Ko(R/Lr) } 
Where Jo and J1 are ordinary Bessel functions and Ko and КІ are modified Bessel functions. 


[Sol] 
a) infinite cylindrical reactor 
i) In the core region 


-DcV"(r) + (г)- уг) 
V^(r) +В2(г)=0, where, В?=(1/К УУ f-Ya5)/D, 
1/r - r- (г) + В?(т)=0 

-(r)=AcJo(BR)+Cc Yo(BR) 


B.C (r)(finite) Cc=0 
(r)=AcJo(BR) 

ii) for the reflector region 

-DcV^(r) + ,(r)=0 

V^(r) -(1/Lr’),(r)=0 


Шт > + <r) - (1/Lr),(r)=0 
,(r)=Arlo(r/Lr)+Crko(r/Lr) 
B.C (r)=0 —Ar=0 
r(r)= CrKo(r/Lr) 
iii) Interface condition 
Фс(К)- Фг(В) 
-Dc <tr) [r=R]=-Dr —(r) [r=R] 


AcJo(BR)+ CrKo(B/Lr)  ................ (1) 
-DcAcBJ,(BR)—drcr91/Lr) K,(R/Lr)..................... 2 
From (2 


Cr=AcJo(BR)/Ko(R/Lr) 

Фс(В)= AcJo(Br) 

Фг(В)= AcJo(Br)/( Ko(R/Lr)): Ko(r/Lr) 

P=K) f 2nrdr 

= 2nk» f Ac-(R/B)-J,(BR) 

Ac-PB/2nk) f В Ј, (ВК) 

(г)= PB/2zk) f R J,(BR)-Jo(Br) 
(г)= PB/2zk) f R J,(BR):-Jo(Br)/Ko(R/Lr)-Ko(r/Lr) 


b) From @, @ 


Jo(Br) Ko(R/Lr) — 0 
-DcBJ,(BR) -Dr(1/Lr)K,(R/Lr) 


Critical condition 
DcB J,(BR)/ Jo(Br)/ = Dr/Lr : K,(R/Lr) / Ko(R/Lr) 


32. The core of an infinite planar thermal reactor consists of a solution of ^"Pu and НО with a plutonium 
concentration of 8.5 g/liter. The core is reflected on both faces by infinitely thick Нь, reflectors. (a) reflector 
savings, (b) critical thickness of the core, (c) critical mass in g/cm’. 


[Sol] 

Z=Ne‘ar/ Nwaa2prMygage 207 C), (Eo) Mg, (Eo) —8.5:107-18-1.0723-1011.3/1-239-0.664 
-1.045 

f2Z/741-0.511 

L4-(1-f) Lry =(1-0.511)8.1=3.96cm? 

1,227 cm? 

Мт^=3.96+27=30.96 cm? 


Reflector savings 
5=7.2+0.1 x (M4?-40.0)26.296cm 


b) 6=0.5(a,-a) — a=a,-26 
Bj^-(u/ a,)’=k -1/ Мт?= 
a,=aV M4j?/ К-1= пу Мт// єтЇ-1=л\30.96/(2.035-0.511-1) 
=87.4cm 
a=87.4-2°6.296=74.8cm 


c) Mp=pr’a= (Ngay Ме/ Nm Мм)а 
= 6.4:10.4(1-239/18) -748 
-0.636g/cm? 


33. The core of a thermal reactor consists of a sphere, 50cm in radius, that contains a homogeneous mixture of 
??U and ordinary H20. This core is surrounded by infinite thick H20 reflector. (a) What is the reflector savings ? 


(b) What is the critical mass? (c) If the maximum thermal flux is 1 x 10P neutrons/cm2-sec, at what power is the 
reactor operating? [ Hint: compute the for core assuming the reactor is bare and of radius 50 cm. Use this 
value of (о estimate 6 and then compute new assuming reactor is bare and of radius 50 + б. Interate 


until convergence is obtained. | 


[Sol] 

a) 

£122.065 , L47-8.1cm?, Ттм-27сш” 

В”-02-( 321-3.948-10?cm? 

7-1- В2(Гл2+ттм)/ £1-1- В “ттм-11-3.948:107(8.1427)/(2.065-1-3.948:107-27)-1.188 
f=Z/Z+1=5.430-10" 

L/’=(1-f) Lyy’=(1-5.430-107)8.1=3.702cm7 

MT?- Lj + ттм = 3.702+27=30.702 cm” 


Reactor savings 
627.240.10(L4? -40.0) 
=7.25+0.10((30-40.0) 
=6.270 


b) 
R=50+ 5=56.270 
В”-(7/56.270)”-3.117 х 10°cm 
-11-3.117:107(8.1-27)/(2.065-1-3.117:107:27)-1.131 
f=Z/741=5.308-10" 
L/’=(1-f) Lyy?=(1-5.308-10)8.1=3.801 
MT?-3.801427-30.801 cm? 
pr7 pw Z-(Mra (Eo) Mygar(20^C) gar(Ko) 
= 1 x 1.131 x 235 x 0.589/(18 x 0.978 х 680.8)-1.306 х 10?g/cm? 
&1=У * pg-4/3-x R3: ру-4/3-л(50) -- ре 
= 6.839kg 


с) 


(г)=А —— . , B=a/R=5.583 x 10?cm'! 

max=(r)=AB 
A=max/B=1-10°/5.583-107=1.791-10'* neutrons/cm:sec 
P=42 AER) ((SinBR-BRcosB R)/B? 


-4л:1.79 x 10 х 3.2:10 x 1.901 x 10?(Sin2.792-2.792c0s2.792)/3.117 x 10° 
-1.303 x 10°W=1.303MW 


34. A spherical-breeder reactor consists of a core of radius R surrounded by a breeding blanket of thickness b. 
The infinite multiplication factor for the core К is greater than unity, whereas that for the blanket k 18 less 
than unity. Using one-group theory, derive the critical conditions and expressions for the flux throughout the 
reactor. 


[Sol] 
At the core region 


-DcV^c +Уас=К Уас 
V’c+B’c=0 where, К Lc?- B? 
c(r)= silBar- —cosB,r 


(r)=finite — A’=0 c(r)= sinB,r 


At the blanket 
-DbV?b -Yabb- К by bb 
V?b«B/?b-0 where B?b--( K b-1/Lb?) 


ь(т)= —вїпВьг+ -eosB,r 


B.C 
p(R+b)=0=,(r)=  sinB,(R+b)+ eesB,(R+b) 
C,=-tanh[B,(R+b)] 
p(r)= —sinByr - —tanh(B,(R+b) cosB,r 
=—sech[B,(R9b) |[sinhB,r cosh B,(R+b)-coshB,r coshB,r 
=—sech B,(R+b):sinh B,(1r-R-b) 


Interface condition 


і) ,(R)=(R) 
A/R-sinBcR=C,/R-sech B,(R+b)sinh(-B,b) 


= -C,/R-sech By(R-b)sinh(B,b)................... (D 


ii) .(r)=A/r?-(B.r-cos B.r-sin B.r) 
г) =С г весв B,(R+b)[ B.rcoshB;(r-R-b)-sinhB,(r-R-b)] 


De ((1)-0, (R) 


D, A[B.RcosB,r-sinB,r]-C;sechB,(R-b)[B.rcoshBy(r-R-b)-sinB;(r-R-b).............. 2 


From the eq. (0,2, the criticality condition is 
sinBcR -sechB,(R+b)sinhB,b 


Dc[BcRcosB,r-sinB,r] sechB,(R+b)[ B.rcoshB,(r- R-b)-sinB,(r-R-b) 


(BcRcotB,b+1)=D,/D,°[1-B. RcotB,R] critical condition 
C1=-A-sinB,R/[sechB,(R+b):sinhB,b] 


-(r)=A/r-sinB,r 
p(r)=- A/r x sinB.R[sechB,(R+b):sinB,b] x sechB,(R+b)sinhB,(r-R-b) 
=- A/r x [sinB, R/sinhB,b] x sinh B,(r-R-b) 


PERY f dv= Еву: B.r-4nr^dr 
= 4 Еу? 
А=р/4 Еву, R? 
c(r)2 p/4 Ев» t В? x sinB,r/r 


ы(г)=- р/4 Eg? + В? x sinB, R/sinhB,b x sinhB,(r-R-b)/r 


35. Plate-type fuel elements for an experimental reactor consist of sandwiches of uranium and aluminum. Each 
sandwich is 7.25 cm wide and 0.16 cm thick. The cladding is aluminum that is 0.050 cm thick. The meat is an 
alloy of fully enriched uranium and aluminum which has 20 w/o uranium and a density of approximately 3.4 
g/cm’. (a) Estimate the mean free path of thermal neutrons in the meat and cladding in this fuel element. (b) Is a 
reactor fueled with these elements quasi-homogeneous or heterogeneous? 


[sol] meat : 20 w/o uranium + 80 w/o aluminum 
p of meat : 3.4 g/cm? 


(a) mean free path à = — 
i) in cladding 
= —= 0 0——— — — — «trem Table П.3> 


= 6.02595 10” atoms/cm? 
6, = 6, + 6, = 0.230 + 1.49 = 1.72 = 1.72 1074 cm’ «from Table II.3» 


he = —=9.64819 =—9.65-em— 
ii) in meat 
N(Al) = 0.8  — — — —«from Table П.3» 
= 6.07081 10” atoms/cm? 
N(U)20.2  »— —— —— ——«frem atom.kaeri.kr> 


1.74222 10? atoms/cm? 
o,(Al) 21.72 1074 cm’ 
o,(U) = 698.2 107 cm? «from atom.kaeri.kr> 


doy = = OE 0960-46 — — — — —— 
(b) à = А. + Am = 10.41 cm 
à >> thickness of a fuel rod(0.61cm) «from 310 page> 
Э quasi-homogeneous 


36. The fuel sandwiches in the reactor in the preceding problem were braised into aluminum holders and placed 
in a uniform array in an ordinary water moderator. The (total) metal-water volume ratio is 0.73, and there are 
120 atoms of aluminum per atom of uranium. Calculate: (a) the thermal utilization; (b) k; (c) the thermal 
diffusion length. [Note: In part (c), compute the value of for the homogeneous mixture by cladding together 


the macroscopic transport cross-sections of the aluminum and water. For 24, of water, use- , where 15 the 


experimental value of for water, corrected of course for density. | 


[sol] (a) Ушеш + Vwater = V 
—=@.73 
V water = 0.578V, V metai = 0.422V 
Манг = — = ------т1934 10° 107 atoms/cm? 


Nmetat = 0.422 — 1.412 10? 102 atoms/cm? 
Nai + № = Nmetat 
—z 120 


Nu = Ngeta = 1.17 10“ 10” atoms/cm? 


Ма = Ма = 1.4 10” 107 atoms/cm? 
f= Lar/ (Zametat + Lawater) = Npoav/ (NuGau + NaiGaal + Nvateroawater) 

= Nugau(T) 6au/(Nugau(D) Gau + Мб! + NwaterGawater) 

= 1.17 10^ 0.978 681/(1.17 10“ 0.978 681+1.4 10? 0.23 + 1.934 10° 

0.664) «from Table 3.2 and Table II.3> 

= 0.829 

(b)k ==: f«Eq.6.77» 
=2.065 0.829 = 1.712 «from Table 6.3> 
(c) D = 44/3 = 1/3X4, «Eq. 5.10» 
Lr = D/X, <Eq. 5.62» 


Y = N waterOtrwater + Моби + М лбу!) 

Ltrwater = 0.58 1/30 = 0.58 1⁄3 0.16) = 2.083 0.58 = 1.208 cm 

Le —41.17 10^ 6.8+1.4 102 3.1) + 1.208 «from Table 6.1> 
-1.25 cm? 

D=1/(3 1.25) = 0.27 cm 

2 = —ЄМ waterÜawater + №620 T NaiGaai) 


= (1.934 10? 0.66441.17 10^ 0.978 681+1.4 10? 0.23) 


«from Table П.3> 

= 0.0833 ст’! 
L4? = 0.27/0.08333 = 3.24 ст? 
Ly; = 1.8 cm 


37. A heterogeneous uranium-water lattice consists of a square array of natural uranium rods 1.50 cm in 
diameter with a pitch of 2.80 cm. Calculate the: (a) radius of the equivalent cell; (b) uranium-water volume ratio 


Vy/Vw; (c) thermal utilization; (d) resonance escape probability; (e) fast fission factor (from Fig. 6.10); (f) k 


[sol] a = d/2 = 0.75 cm 
(a) tb” = 2.807 
b = 1.58 cm 


(b) = ——=0.291 


(c) = —+Е 
Е(х) = — + - 2. —*® —5........ 


Е(у, 1) = —— 


Lr = 1.55 cm, Ly = 2.85 cm 
x = a/Ly = 0.75/1.55 = 0.484 
y = a/Ly = 0.75/2.85 = 0.263 
Z = ЪЛ м = 1.58/2.85 = 0.554 





F(x)=1+-—*- — —*+ — —$ 
= 1.029 

Е(у, 2) =1+ Е - -+ > 
= 1.0411 


Уем = Мм _6в.м(Ео)( +” 
Lar = Ny-Bar(T)oam(Eo)( 2^ 
Nu = ———= №( —— 


-2.59 10? 10^" atoms/cm? 
№= ———= Nx — 

= 1.089 10? 10” atoms/cm? 
№ = 0.00723 = 7.84 10° 10” atoms/cm? 
М» = 0.99283 = 1.081 10? 10^ atoms/cm? 


At 20°C, g,;s( T) = 0.9780, g428(T) = 1.0017 
бам = 0.664 b, 6a25 = 681 b, саз = 7.59 b 


= FTE 
= 1.496 
f = 0.669 
(d) P = exp[-——— 
А = 2.8, C = 38.3, өм? м 71.46, а = 0.75 cm, р = 19.1 g/cm3 
1=А+С/ = 2В + 38.3/ =12.92 
Р = 0.9723 
(е) £ = 1.2 «From Fig. 6.10> 


(f К =efPe= 2.065 0.669 0.9723 1.2 
= 1.612 


38. Repeat the calculations of Problem 6.37 for the case in which the uranium is enriched to 2.5 w/o in U. 


[sol] (a) b = 1.58 cm 

(b) — 3.44 

(c) x = 0.484, y = 0.263, z = 0.554 
Е = 1.029, Е = 1.0411 
№ = 2.59 10° 10” atoms/cm? 
№ = 1.089 10? 10" atoms/cm? 
а/о = [(2.5/235 + 97.5/238) 0.01]°1/235 = 2.5 
№ = 0.0253 = 2.7 10“ 10” atoms/cm? 
Nog = 0.9953 = 1.058 107 10%” atoms/cm? 


= 04 н — — — — — — —— 


= 1275 
f = 0.784 
(d) P= exp[-4= 9.9723 
(e) € =1.2 
(РК “єїїРЕ-2.065 0.784 0.9723 1.2 
= 1.36 


39. In a hexagonal lattice (also called a triangular lattice), each fuel rod is surrounded by six nearest neighbor 
rods, equally spaced at a distance equal to the pitch s of the lattice. Show that the radius of the equivalent cell in 
this lattice is given by 


b=0.525 s. 


[sol] 


с = 60° 
ль? =2 s’sino 
‚ - 
= 98$ — 


b = 0.525 $ 


40. Calculate К for a hexagonal lattice of 1.4 cm radius natural uranium rods and graphite if the lattice pitch 18 
20 cm. [Note: The fast fission factor for this lattice 1s 1.03.] 


[sol] Hexagonal lattice = natural uranium + graphite 
а = 1.4 ст, s = 20 cm, = = 1.03, єт = 2.065 «from Table 6.3» 


i) = (> амУм / У ar Vp)F(x) + E(y,Z) 
x=a/Lp=1.4/1.55 = 0.903 
у= а/м = 1.4 / 55 = 0.0255 
Z=b/Ly = 0.5253 / Ly = 10.5 / 55 = 0.191 
F(x) = xla) / 2 (Хх) 


1 (0.903) = 1.22 

1, (0.903) = 0.498 

Е (0.903) = (0.903 1.22)/(2 0.498) = 1.106 
1 (0.0255) = 1.0006, 1,(0.0255) = 0.013 


1. (0.191) = 0.1, Ko(0.0255) = 10 

К. (0.0255) = 40, К, (0.191) = 4-9 

Е(у, 7) = (22 - у”) / (2у)(16(у) Ki(z) + Ko(y) 112) ES C2 Ki(y) - Ki] (у) 
E(0.0255, 0.191) = 1.654 

Yam / Угар = №с.м(Е0) / Nroar(T)o,r(E0) 

№ = РмМАУм / МмуУ = 3.286 107 10” atoms/cm? 

№ = PrNaVe/ МУ = 8.592 10“ 107 atoms/cm? 

№: = 0.00723 26.186 10° 102 atoms/cm? 


М» = 0.99283 = 8.53 10“ 107" atoms/cm? 
$a25 = 0.978, 9128 = 1.0017 


О.м(Ео) = 0.0034 b, OP; ( E4) = 681 b, Os ( Eq) =7.59b 
dam / Žar Ммбам(ЕО) / Мев.к(Т)ак(Е0) 
= 0.0105 
Ук / Ум = a? / (b? - а?) = 0.0181 
= (0.0105 / 0.0181) 1.106 + 1.054 = 1.696 
f = 0.59 
ii) А = 2.8, С = 38.3, а = 1.4, p = 19.1, этУсм = 0.0608 


I=A+C/ = 2.8 + 38.3/  -10207 


Р = exp[-1-2-0:997 
К = атр = 2.065 0.59 0.997 1.03 
= 1.25 


Chapter 7 


The Time-Dependent 
Reactor 





1. Compute the prompt neutron lifetime for an infinite critical thermal reactor consisting of a homogeneous 
mixture of ^?U and (a)D2O, (b) Be, (c) graphite. 


[sol] Assume room temperature T=20°C 
Table 6.3 


Table 7.1 


(a) 
(b) Be 


(c) Graphite 


2. Calculate the thermal neutron diffusion time in water having a density of 


[sol] Table II.3 


3. Express the following reactivities of a 2?Pu-fueled thermal reactor in dollars (a)0.001, (b) 496, (c) -0.01 


[sol] Table 7.2 


(c) = 
4. Express the following reactivities of a 253U. fueled thermal reactor in percent: (a) 0.001, (b) $2, (c) -5Ocents 


[sol] Table 7.2 


(a) 0.001=0.1% 
(b) - 
(c) -50cents=-0.5dollars= 


5. Plot the reactivity equation for one-delayed-neutron group for а ^?U-fueled thermal reactor and prompt 
neutron lifetimes of (a) Osec, (b) 10^sec, (c) 10?sec. Take 


[sol] | — —— — —— 


(a) 


(b) 


(c) 


6. From the plot in Problem 7.5, determine the periods of ^?U-fueled thermal reactors with ара 
reactivities of (a) 40.196, (b) -10cents, (c) 44.50, (d)+$1.00. 


[sol] 


(a) 


(b) —— 


(c) — 


(d) — 


7. Compare the answers to Problem 7.6, which are based on the one-delayed-group model, with periods 
determined from Fig.7.2. 


[sol] (a) 40.146 , 
: Fig. 7.2 - - 


(b) -10cents 
Fig.7.2 - - 


(c)$0.50 


(d) $1.00 


8. А ^?U-fueled reactor originally operating at a constant power of 1 milliwatt is placed on a positive 10-minute 
period. At what time wil the reactor power level reach 1 megawatt? 
[sol] T=10minutes=10 =600s 


9 


9. Calculate the ratio of the concentration of the 22.7-sec precursor to the thermal neutron density in a critical 
infinite thermal reactor consisting of a mixture of 235U and ordinary water, given that the thermal flux is 1013 
neutrins/cm2-sec.[Hint: place dC/dt =0 in eq. 7.21] 


[Sol] 
dc/dtZ[pK5 a4j/P]-AC  , С=ВКУат/РА, 


р-0.0065, К=єтї=1‚ 4-[o/Vz]-z-v1,v1-2200m/sec 


{= / = / + — -[f/1 -f]- 
{=1/<т=1/2.065=0.484 

=[f/1-f]- ==[0.484/1-0.484]-0.0197=1.85 x 10°?сп! 
Уа-1.85 x 107-0.0197-3.82 х 10? cm” 
A=In2/22.7=3.054 x 10?sec'! Р-1 


C=(0.0065 х 1 х 3.82 x 10 “х (а/п) £ x 2200 х 107/(1 x 0.054 х 10?) 


C/ є-2.018 х 10° 


10. The first term on the right hand side of Eq.(7.19) gives the number of prompt neutrons slowing down рег 
cm3/sec while the second term gives this number for the delayed neutrons. Compare the magnitude of these two 
terms in a critical reactor. 


[Sol] 
First term  (1-B) КУат 
Second term PAC 


From eq. 7.20 
dC/dt- [ВКУат/Р]- С 


1С- РК»аүР 
(1-8) KY.ay-PAC=(rB) КУат-ВКУат 
=(1-2B) ВКУат term Second term 


11. Fifty cents in reactivity is suddenly introduced into a critical fast reactor fueled with 235U. What 18 the 
period of the reactor? 


[Sol] 
p=50cents=0.5 x 0.0064=3.2 x 107 
T=0.0848/3.2 x 1035ес 


12. The reactor in problem 7.8 is scrammed by the instantaneous insertion of 5 dollars of negative reactivity 
after having reached a constant power level of 1 megawatt. Approximately how long does it take the power level 
to drop to 1 milliwatt? 


[Sol] 
Ро-1 x 10^ watt, P=1 х 10^watt 
p=-5$=-0.0325 


Prompt drop P,=Po-B(1-p)/(B-p)=[1 х 10° x B(1+5 х 0.0065)]/(p--5p)—0.1721 Mwatt 
P= P; 

1 x 10?20.1721 x 10* 

T= -80In(1 x 10?)/( 0.1721 x 10°)=1517sec 


13. When a certain research reactor operating at a constant power of 2.7 megawatts is scrammed, it is observed 
that the power drops to a level of 1 watt in 15 minutes. How much reactivity was inserted when the reactor was 
scarmmed? 


[Sol] 

Р,-Р, 

Р-Р,-1 x exp(15 x 60/80) -7.688 х 10^Watt 

Ру= P, -B(1-p)/(B-p) 

P,(B-p)= B(1-p) Р, 

p (P,- Р„)= p(P,- B Pj), p= p (Р,- P,)/ Pi- p P,=0.0065(7.688 x 10*-2.7 x 10°) 
--0.2874 


14. An infinite reactor consists of a homogeneous mixture of 235U and H20. The fuel concentration is 5% 
smaller than that required for criticality. What is the reactivity of the system? 


[Sol] 

From the problem(7-9) 

У`ав=0.0185ст? 

Fuel concentration is 5 % smaller than criticality 


Yar=0.0185 x 0.95=0.0176ст? 


f=) a Y art ar=0.0176/(0.0176+ 9.0222)=0.472 


K-nf=2.065 x 0.472=0.975 
p= K / K=0.975-1/0.975=-0.0256=-2.56percent. 


15. One cent worth of either positive or negative reactivity places a235U-fueled reactor on what period? 


[Sol] 
р=1сепіѕ=0.01 x 0.0065=6.5 х 107 


T=1/p x У;В:6=0.0848/6.5 x 10°sec=1304.6sec 


16. It is desired to double the power level of a 235U-fueled reactor in 20 minutes. If Ip = 10~ sec, (a) on what 
stable period should the reactor be placed? (b) How much reactivity should be introduced? 


[Sol] 
2Po=Po, 1,=10“sec 
a) 20/T=In2 
T=20/In2=28.85 min=1731sec 


b)Using Fig.7.2 

1,=0.0001sec 

T=1731sec 
p—cannot be found very small reactivity 
Using Eq7.40 


T-0,7) iiti) 
p= d,tY;pit)- 440540.0848)-4.905 x 10?—0.007546p 
0.7546cents 


17. During test-out procedures, a 239Pu-fueled thermal reactor is operated for a time at a power of 1 megawatt. 
The power is then to be increased to 1 00 megawatts in 8 hours. (a) On what stable period should the reactor be 
placed? (b) What reactivity insertion is required? 


[sol] 
(a) 


T = t/In(P/ ——— 


(b) 


19. An experimental reactor facility is a bare square cylinder 100 cm high, composed of small beryllium blocks 
with thin foils of 235U placed in between, so that the system can be considered to be a homogeneous mixture of 
Be and 235U. The reactor is to be controlled with a single black control rod 2.5 cm in radius and located along 
the axis of the system. (a) If the reactor is critical with the rod fully withdrawn, how much negative reactivity is 
introduced into the system when the rod 18 fully inserted? (b) Assuming that the rod moves into the reactor 
instantaneously, on what period does the reactor go? 


[sol] 
(a) 


d = 2.131x0.5x 


where 


— —=8.1112 


T= ----- —33х зе _— 


20. Suppose it is desired to control the reactor described in Example 7.7 with one central 
rod having a worth of 1 0%. How big should the rod be? [Hint: Plot Pw versus a.] 


[sol] 





= 0.163[0.116--In( -= 0.1 


In( -= 1414 





d- 2.131x0.16x 


In —9-4 =b444— 


a=4.3 cm 


2] . If the reactor 1 n Problem 7. 1 9 i s controlled b y an array o f 2 5 more or less uniformly distributed black 
rods 0.3 in. in diameter, what is the total worth of the rods? 


[sol] 
а = 0.381cm 


d = 2.413cm 


y = 0.381/98 = 0.0389 


7. = 10/98 = 1.03 


E(y,z) шал 


= —= 6.785 => = 0.147 


22. A certain pressurized-water reactor is to be controlled by 61 cluster control assemblies, each assembly 
containing 20 black rods 1 . 1 5 cm in diameter. The reactor core is a cylinder 320 cm in diameter. The average 
thermal diffusion length in the core is 1 .38cm, D = 0.2 1 cm, and Lt in the core material is approximately 2.6 
cm- I Calculate the total worth of the rods. 


[sol] 


а= — 


уза/-0417, Z= =3.320 
E(y,z) = 2.625 





=> — => —4-0:5425 + 2.625 = 7.742 = 0.1292 


23. Consider the problem of detennining the worth of а cruciform control rod in the onedimensional geometry 
shown in Fig. 7. 1 0. First suppose that qT fission neutrons slow down to thermal energies per cm3/sec. The 
quantity qT can be treated as a source term in the thermal diffusion equation, 

(a) Show that the solution to this equation, subject to the boundary conditions, 

(b) Compute JR by multiplying the neutron current density at the blade surface by the area of the blades in the 
cell and then dividing by the total number of neutrons thennalizing per second in the cell 


[sol] 
(a) 
D— 


(b) 
J(m/2 - a) = -D — a ee 


24. It is proposed to use the crucifonn rods described in Example 7.9 to control the PWR described in Problem 
7.22. How far apart should these rods be placed to provide the same total worth as the cluster control rods in that 
problem? 


[sol] 
В = 320/2 = 160 cm, =1.38cm D=0.21cm, ‚ »D=0.4475 =0.1723 


= — —=9.147, а = 0.312 х 2.54 2=0.396 cm, 129.75 x 2.54/2 = 12.38 ет 


0.147 = ———— к (0.3244coth —)=-450 


=> т = 192cm 


25. Suppose the fast reactor described in Example 6.3 is controlled with 50 rods, each rod containing 
approximately 500 g of natural boron. Estimate the total worth of the rods. 


[sol] 


26. The core of a fast reactor 18 a square cylinder 77.5 cm in diameter. The composition of the core by volume 18 
as follows; 25% fuel; 2596 stainless steel cladding and structural material; 50% liquid sodium. The fuel consists 
of a mixture of 238U and 239Pu having a density of 193.1 g/cm3, the plutonium making up 20 w/o of mixture. 
It is desired to provide 5.5% reactivity control in 10 control rods containing В4С. How much ВАС 1$ required 
per rod? 


[sol] 


27. A control rod 100cm long is worth 50 cents when totally inserted. (a) How much reactivity is introduced into 
the reactor when the rod is pulled one-quarter of the way out? (b) At what rate is reactivity introduced at this 
point per cm motion of the rod? 


[sol] 


28. Suppose that, at some time during its operating history, the reactor described in Example 7.7 is critical with 
the rod withdrawn one-half of its full length. If the rod is now suddenly withdrawn another 10 cm. (a) hoe much 
reactivity is introduced? (b) on what period does the reactor power rise? 


[sol] 


29. What concentration of boron in ppm will give the same worth as the single control rod in Example 7.7? 


[sol] 


30. An infinite 235U-fueled, water-moderated reactor contains 2096 more 235U than required to become critical. 
What concentration of (a)boron in ppm or (b) boric acid in g/liter is required to hold down the excess reactivity 
of the system? 


[sol] 


3]. A reactor is operating at constant power. Suddenly the temperature of the incoming coolant drops below its 
previous value. Discuss qualitatively the subsequent behavior of the reactor in the two cases; (a) aT is positive, 
(b) aT is negative. Show, in particular, that the reactor is unstable in (a) and stable in (b). 


[sol] 
017| 0| 41 z 9j оны A[9 
32. Calculate the prompt temperature coefficient at room temperature of a reactor lattice consisting of as 


assembly of 1-in diameter natural uranium rods in a heavy-water moderator, in which the moderator volume to 
fuel volume ration is equal to 30. 


[sol] 


33. А C0»-cooled, graphite-moderated reactor fueled with rods of slightly enriched uranium dioxide 1.5 cm in 
diameter has a resonance escape probability of 0.9 1 2 at 300K. What is the value of p at the fuel operating 


temperature of 665°C? 


[sol] — 


34. The overall temperature coefficient of a ^?U-fueled reactor is -2x 107 per °C and is independent of 
temperature. By how much does the reactivity of the system drop when its temperature is increased from room 


temperature (about 70°F) to the operating temperature of 550°F? Give your answer in percent and in dollars. 
[Note: The decrease in reactivity calculated in this problem is called the temperature defect of the reactor.] 


[sol] — m == 


_ 20.994684 


35. Show that the temperature coefficient of B^ is given by 


where [В is the coefficient o f volume expansion o f the reactor structure. [Hint: Recall that В is equal to three 
times the coefficient of linear expansion.] 


[sol] —— 


1) For spherical reactor 


36. A pressurized water reactor fueled with stainless steel-clad fuel elements is contained in a stainless steel 


vessel that is a cylinder 6 ft in diameter and 8 ft high. Water having an average temperature of 575^F occupies 
approximately one-half of the reactor volume. How much water is expelled from the reactor vessel if the 


average temperature of the system is increased by 10°F? [Note: The volume coefficients of expansion of water 


and stainless steel at 300°C are 3х 107 рег °C and 4.5x 10? per °C, respectively.] 


[sol] - 


37. The thermal utilization of the reactor in Problem 7.36 is 0.682 at 575^F. Using the results of that problem 


and ignoring the presence of structural material in the core, estimate O7(f) at approximately 575°F. 


[sol] — 


38. What is the effective half-life of Xe in a thermal flux of 10'* neutrons/cm?-sec at a temperature of 800°C? 


[sol] — 


39. Compute and plot the equilibrium xenon reactivity as a function of thermal flux from 
Q1 = 5x10” to Фт = 5x10! 


[sol] ——— 


For u-235, 








1Е+12 1E+13 1Е-14 1Е+15 

















40. Using Fig. 7. 14, plot the maximum post-shutdown xenon reactivity as а function of thermal flux from 
Фт = 5x10! to Фт = 5x10! 


[sol] 





1Е}12 1Е+13 1E+14 1Е+15 


41. A ^?U.fueled reactor operating at a thermal flux оЁ5 10 neutrons/cm-sec is scrammed at a time when 
the reactor has 5% in reverse reactivity. Compute the time to the onset of the deadtime and its duration. 


[sol] 
- d T + Ух)Фт e^ a Vir (e _ о-и) 
vp Фх + ÖT Pi — Px 
<from Eq. 7.103> 
у = 2.42 
pe=1 


үт + yx = 0.066 «from Table 7.5» 

y; = 0.0639 «from Table 7.5» 

Өт=5 10" 

Øx 20.770 10” «from Eq. 7.97» 

©, = 1.055 107 «from Eq. 7.104» 

dx = 0.0753 hr-1 «from Table 7.6» 

№ = 0.1035 hr-1 «from Table 7.6» 

р = -0.486877e °°7°°* + 9.463244 1055 
= -0.05 

t 2 3.20253 and 19.7751 

the onset of the deadtime - 3.20253 hr 

duration = 16.5726 hr 


42. Calculate the equilibrium concentration, in atoms/cms, of Уе and '?Sm in an infinite critical ^"U-fueled, 
water-moderated thermal reactor operating at a temperature of 200°C and a thermal flux of 10" neutrons/cm’- 
sec. 


[sol] 
43. How much reactivity is tied up in '?Xe and "Sm in the reactor described in Problem 7.42? 
[sol] 


44. Gadolinium-157 is a stable nuclide having an absorption cross-section at 0.0253 eV of 240,000 b. It is 
formed from the decay of the fission product P'Sm according to the following chain: 


157% т > > 
1 


Fission 


Neither "Sm nor Eu absorbs neutrons to a significant extent. The ^ `U fission yield of P'$m is 7 10? 
atoms per fission. (a) What is the equilibrium reactivity tied up in "Gd in a reactor having an average thermal 
flux of 2.5 10" neutrons/cm^-sec? (b) What is the maximum reactivity due to this nuclide after the shutdown 
of the reactor in part (a)? 


[sol] (a) — r,Y От - XS 9 S = гу, Or/ А, 

== AS " ЕЕ Э Е = AS / kk 51.1 От / Ар 

— ЕЕ - o4 ;00G 
Э G = Е /с„Өт = г.у, От / в, От = г / бас 

р= - (У.с / УВ / Ур== - (С сс; / У) /урғ = - г, /урё=-(7 10°) / 2.418 
= - 2.895 10° 
=- 2.895 10° % 

(b) G(t) =С,+Е А е^) 
If Саара E, are equilibrium concentraties prior to shutdown 
Go = G , E, = Е 
p(t) = - Qaa / >) / VPE = -(басС(0/ 24) | ур = - (1 / vps)(r, + (г,басӨт / №) - ё? 2) 
= - (г, / vpg)(1 + (Ør/ Og)(1 - е^") 
Р ах = == (r, / урғ)(1 + (Әт/ 06) 
Ок=к/ бас-( /(1.52 3600)/(2 240000 1075 
-5.955 10" atoms/cm?sec 


45. A nonradioactive fission produce has an absorption cross-section of 75 b. Should this be considered a 
permanent poison in a reactor having a thermal flux of 3 10" neutrons/cm"-sec? 


[sol] 


10° t= 105 neutrons/cm?-sec 


l 3х 104 - 


10/4 









Poison cross-section, barns/atom 235 U 


1012 


1 anul 
105 10° 107 10° 


Irradiation time, sec 


Figure 7.18 Buildup of fission product poisons in bams of absorp- 
tion cross-section per original atom of 2350. (From ANL-5800, 2nd ed., 
1963.) 


t= 62 10's 


This is considered a permanent poison. 


46. What fraction of the poisoning in Example 7.14 1s due to Xe? 


[sol] p. = - —Өт/(Өх + 01) = - 
= -3.156 10° 


in Ox = Ax / бах = = 
= 7.684 10" atoms/cm?sec 
From Example 7.14 
= -0.00849 


fraction of the poisoning = — ——=37-47 % 


47. An infinite reactor containing no fertile material operates at constant power over its lifetime. (a) Show that 
the atom density of the fuel decreases according to the relation 


М (0 = N&(O)[1 -„„@т(0){], 


Where №(0) and @т(0) are, respectively, the fuel atom density and thermal flux at startup. (b) Find an 
expression for the flux as a function of time. 


[sol] (a) №(0) = №0) - O4(0)dt 
= №(0) — Np(Q}O7(0)t 
= №(0) [1 ғӘт(0) 0 
(b) Yar(t)Or(t) = У,к(0)01(0) = Constant 
Dar(t)Or(t) =ғ№(0)(1 -„к@т(0)] Or(t) 
= $ar(0)[1 -arØr(0)t] Or(t) 


= Žar(0)Ør(0) 
От(0 = O4(0) / (1 -кОт(0)0 


48. An infinite 139Pu-fueled fast breeder has the breeding gain С and 18 operated at constant power before 
refueling. Derive expressions for the concentration of 239Pu and the (one-group) fast flux as a function of time 
after startup. 


[sol] Pu?" + 5 fission and Pu?" 
Concentration of "Pu 
= -№с О(г, t) — №с:Ө(г, t) 
= -N»«(6; + 6,)O(r, t) 
where, o; = fission cross section 
c, = radioactive capture cross section 








Chapter 9 


Heat Removal from 
Nuclear Reactors 


]. The nuclear ship Savannah was powered by a PWR that operated at a pressure of 1,750 psia. The coolant 


water entered the reactor vessel at a temperature of 497°F, exited at 519°F, and passed through the vessel at a 


rate of 9.4 x 1 06 lb/hr. What was the thermal power output of this reactor? 


[sol] 


2. An experimental LMFBR operates at 750 MW. Sodium enters the core at 400°C and leaves at 560°C. At what 
rate must the sodium be pumped through the core? 


[sol] 


3. Show that if water passes through the core of a BWR at the rate of w lb/hr, then уу, lb/hr of steam will be 
produced, where 


q is the thermal power of the reactor; h; and hin are, respectively, the specific enthalpies of the saturated water 
leaving the core and the water entering the core; and һе, is the heat of vaporization per pound. 


[sol] 


4. ABWR operates at a thermal power of 1,593 MW. Water enters the bottom of the core at 526°F and passes 
through the core at a rate of 48 x 1 06 Ib/hr. The reactor pressure is 1,035 psia. Using the result of the previous 
problem, compute the rate in Ib/hr at which steam is produced in this reactor. 


[sol] 


5. Starting with the value of the specific enthalpy of saturated water at 350°F given in Table IV.1, calculate the 


specific enthalpy of saturated water at 600°F. [Hint: Integrate the specific heat using Simpson's rule.] 


[sol] 


6. The decommissioned Fort St. Vrain HTGR was designed to produce 330 MW of electricity at an efficiency of 
39.23%. Helium at a pressure of 710 psia enters the core of the reactor at 760°F and exits at 1,430?F. What is the 


helium flow rate through the core? 


[sol] —— 


7. А small PWR plant operates at a power of 485 MWt. The core, which is approximately 75.4 in. in diameter 
and 91.9 in high, consists of a square lattice of 23,142 fuel tubes of thickness 0.02 1 in and inner diameter of 
0.298 in on a 0.422-in pitch. The tubes are filled with 3.40 w/o-enriched UO;. The core is cooled by water, 


which enters at the bottom at 496°F and passes through the core at a rate of 34x 10° Ib/hr at 2,015 psia. Compute 


(a) the average temperature of the water leaving the core; (b) the average power density in kW/liter; (c) the 
maximum heat production rate, assuming the reactor core is bare. 


[sol] a) 


9. The core of a BWR consists of 764 fuel assemblies, each containing a square array of 49 fuel rods on a 0.738- 
in pitch. The fuel rods are 175 in long, but contain fuel over only 144 in of their length. The outside diameter of 


the fuel rods is 0.563 in, the cladding (fuel tube) is 0.032 in thick, and the ОО, fuel pellets are 0.487 in diameter. 


This leaves a gap of (0.563 - 2 x 0.032 - 0.487) / 2 = 0.006 in between the pellets and the cladding. The UO; 
has an average density of approximately 10.3 g/cm’. The radius of the core is 93.6 in, and the reactor is designed 
to operate at 3,293 MW. The peak-to-average power density is 2.62. Calculate for this reactor: (a) the total 
weight of 002 and uranium in the core; (b) the specific power in kWikg U; (c) the average power density in 
kW/liter; (d) the average linear rod power а ay in kW/ft; (e) the maximum heat production rate. 


[sol] 


a) = 


9. The core of an LMFBR consists of a square lattice of 13,104 fuel rods 0.158 in in diameter, 30.5 in long, on a 
0.210-in pitch. The fuel rods are 26 w/o-enriched uranium clad in 0.005-in stainless steel. Liquid sodium enters 
the core at approximately 300°C and passes through the core at an average speed of 31.2 ft/sec. The core 
produces 270 MW of thermal power, with a maximum-to-average power density of 1.79. Calculate: (a) the 


maximum heat production rate; (b) the maximum neutron flux. 


[sol] (а) ч = 1.16PE, / Ha?nEg 
= (1.16 270 10° 3412« . ——-— 180<MeV>) 


/[30.5«in» (0.158/2)4 40» 0.000579<  »— 13104 200<MeV>] 
= 6.663 10° Btu/hr-ft 
(b) p(U) = 19.1 g/cm? 
№0025) = (0.26 19.1 0.6022 102)/235-0.0127 10^ atoms/cm? 
Уһ = 0.0127 582.2 = 7.41 cm” 
-(/2) 0.931 7.41 ( 32 24:37 cm" 
Ømax = O(r=0, 7-0) = q шах/( Еа) 
= (6.663 10°<Btu/hr-ft?> 6.461 10’< —>)/4.37<ет"> 180<MeV>) 
= 5.47 10" neutrons/cm?-sec 


10. The variation of the neutron flux and/or heat production rate along the 7 direction in the core of a reflected 
reactor is often approximated by the function 


Ø = constant cos( >, 


where is the distance between the extrapolated boundaries and is somewhat larger than Н, the actual height 
of the core. Show that the maximum-to-average heat production ratio in the z direction Q, is given by 


Q, = 


[Note: A similar approximation can be made for @ and/or 9" in the radial direction by writing 0 = constant 
Jo(2.405r / ).] 


[sol] лах = О(7:0) -С 
О, = / = — — 
= -sin( > 
О, = дһ. / Day = C/ -sin( — = 


11. Show that the maximum-to-average flux and/or power distribution ratios in the axial and radial direction are 
related to the overall ratio by 


Where the notation 18 that of Problem 8.10. 


[sol] Ø = CJo( —5 
О, = ээрээ 
In general, 
Or, z) = CJo—93cos( > 


Øn = - --- 2 
О ах = CJo(0)cos(0) =C 
Q = О, ах / О., - 


= О, О, 


12. The соге of a fast reactor is а cylinder 38.8 cm in radius and 77.5 cm high. Two dimensional (т, 7) 
multigroup calculations show that the power density distribution in the core can be represented approximately 
by the expression 


P(r, 2) = Po1— 0-]соз( 3— 


Where Po is a constant and r and z are the distances in centimeters from the axis and the midplane of the core, 
respectively. (a) Evaluate Ро in terms of the total core power Р. (b) What is the maximum-to-average power ratio 
in the core? (c) Calculate the maximum-to-average power ratios in the radial and axial directions. 


[sol] (a) Р- 

= 2Р, = — шин 
- Р, —— - 
= PR'- 65 — = 
= 109P, [38.87 - € 7] 2sin( + 
=2.096 10°P, 

Po = 433+ 10°P 

(b) Pos = Р(т, Z)r=0, z=0 = Po 
Pus = S 


= 0.572P, 
О=—= —-=k748 


(c) P(r) = P4[1 — 07] 
P(r)max = Pro 
P(r), = —— 


=Р[1- € X] 
= 0.7106P,, 
О. = P(r) max / Pro = —4:407 


OQ = OQ, О, 
О, =О/О, = 1.748 / 1.407 = 1.242 


13. ABWR power plant operating at an efficiency of 34% has an electrical output of 1,011 MW. (a) What is the 
maximum fission product decay energy in the reactor at shutdown? (b) What is the decay energy 6 months after 
shutdown? 


[sol] (a) Po = 1011 / 0.34 = 2973.53 MW 
At shutdown, the shortest time in Fig. 8.3 is 107 sec and the longest time in Fig. 8.3 is 10° sec. 








0.07-3 10“ 
= 0.0697 
Рак = 0.0697 Р,-207.255 МУУ 


10 10% 10! 102 10? 104 





Infinite reactor operating time 








107 10° 10° 10? 108 10° 


Time after shutdown (cooling time), !,, sec 


Figure 8.3 The ratio P(1,)/ Ро of the fission product decay power to 
the reactor operating power as a function of time г, after shutdown. (Sub- 
committee ANS-5 of the American Nuclear Society, 1968.) 


(b) t, = 6 months = 1.58 107 sec 
If the reactor operates for one year, tọ = 3.16 10’ sec 


— 
— 


= 0.001 — 0.0007 = 0.0003 
P=0.0003 2973.53 = 0.892 MW 


14. А 1,000-MW reactor is operated at full power for 1 year, then at 10% power for 1 month, and the shut down. 
Calculate the fission product decay heat at shutdown and for 1 month after shutdown. 


[sol] (a) at shutdown 
(i) For 1,000 MW 
t, = 1 month = 2.6 10° sec 
to = 1 year = 3.16 10’ sec 


———Á 
— 


-18 10°-6.8 10“=1.12 10° «From Figure 8.3» 
P=1,000 1.12 10°=1.12MW 
(ii) For 100 MW 
t, = 107 sec 
to = 1 month = 2.6 10° sec 


=0.07-1.8 10° = 0.0682 «From Figure 8.3» 

P=100 0.0682 = 6.82 MW 

the fission product decay heat at shutdown 

= 1.12 MW + 6.82 MW = 7.94 MW 

(b) for 1 month after shutdown 
(i) For 1,000 MW 

t,=2 month = 5.2 10° sec 
to = 1 year = 3.16 10’ sec 





— 


-15 10°-7.2 10°=7.8 10“ <From Figure 8.3» 
P=1,000 7.8 10“=0.78 MW 
(ii) For 100 MW 
t,=1month=2.6 10° sec 
to = 1 month = 2.6 10° sec 


-18 10?-1.5 10°=3 10^«From Figure 8.3» 
P=100 3 10*20.00MW 
the fission product decay heat at shutdown 
= 0.78 MW + 0.02 MW = 0.8 MW 


15. A power reactor fueled with ^?U operates at the power Р, for 1 year and is then scrammed by the insertion 
of $9.63 in negative reactivity. (a) Calculate the level to which the fission power of the reactor immediately 


drops. (b) To what power level does the reactor actually drop? (c) Compute and plot the reactor power from 0.1 
sec to 0.5 hr after shutdown. 


[sol] (a) p =$9.63=9.63 0.0065 = 0.06259 
The fission power drops suddenly t the level P4, where 
Р, = —Ру- 
= —P,--9:09997p, 
(b) The fission produce decay heat 





—É— 





= 0.07 — 0.0007 = 0.0693 «From Figure 8.3» 
P, = 0.0693P0 


Since the actual power are produced by the fission within the reactor and the decay of radioactive 
fission products, the actual power is given by 

P = 0.09997P, + 0.0693P = 0.16927P, 

By fission 

Р; = Pye" = 0.09997P,e 5° 

Where T = 80 sec 

By the decay of fission product 


Ра=— ---Р, 
Р = 0.09997Р,е", 89 +(—— ---Р, 


(i) t, = 0.1 sec 
P = 0.09997P,e *5* (0.07 —0.0007)Р‹ 
= 0.16915P, 


(ii) t, = 0.5 hr = 1,800 sec 
P -0.09997Р,е! 099,---- -----Р, 


=1.69 10°P,+(1.8 10°- 0.0007)P, 
-0.0173Р, 


(c 


Хэ 





16. A thermal reactor fueled with slightly enriched ^"?U is operated at a power of 2,800 MW for a period of 2 
years and then shut down for refueling. The reactor has a conversion factor of 0.82. (a) What 18 the fission 
product decay energy at shutdown and 1 month later? (b) How much decay heat is due to ^"U and ^"Np? (c) 
What is the activity of the fission products and the ^"U and ^PNp at the prior times? 


[sol] (a) (i) at shutdown 


-0.07-6.5 10“ «From Figure 8.3» 
= 0.06935 
P-2800 0.06935 = 194.18 MW 
(ii) one month later 


—— m— 


=1.8 10“-6.45 10“ <Егот Figure 8.3» 
= 1.155 10° 
P=2800 1.155 10°=3.234MW 
(b) (i) For 2710 
a. at shutdown t, = 0 
—=2.28 10° ( ЭН-| 


=2.28 10° 0.82 ( = 
=2.187 10° 

Р, = 2800 2.187 10° = 6.12 МУУ 
b. one month later 


Р» = 0 
(ii) For ^?Np 
a. at shutdown t, = 0 
—= 2.17 10? ( HA - 1 


-7.0 10°(1- ) } 
=2.17 10° 0.82 ( 3— 
= 2.082 10° 


Р, = 2800 2.082 10? = 5.83 MW 
b. one month later 


— 2.082 10? =2.94 107 


Р, = 2800 2.94 107=8.3 10“ МУУ 
(c) (1) fission product activity 

-14 10*P[t;'? — (t, +) ^] 

a. at shutdown, t, = 0.1 sec = 1.16 10° day 
fission product activity 
=1.4 10°P[(1.16 109) (2 365)'?] 
=2.11 10’ 2800 
=5.91 10” 

b. one month later, t, =30 day 
fission product activity 
=1.4 10°P[(30)°?-(30+2 365)9 | 
= 3.38 10° 2800 
=9.45 10° 

(ii) For ^?U 

a. t, = 0.1 sec 
fission product activity 
=5.91 10" —-k86 10° 


b. t, = 30 day 
fission product activity 
=9.45 10° —2.98 10’ 


(iii) For ??Мр 


a. t, = 0.1 sec 
fission product activity 


=5.91 10" --577 10° 


b. t, = 30 day 
fission product activity 


=9.45 10° —=284 10’ 
17. The reactor core in Problem 8. | 2 produces 400 MW, of which 8 MW is due to y-ray heating of the coolant. 
The total heat transfer area of the fuel is 1,580 ft2 Compute (a) the average power density in the core in kWniter 
and kW/ft3 ; (b) the average heat flux; (c) the maximum heat flux. 


[sol] 


(a) 


= 1.07mw/l = 1070kw/l, where P(1,z) => [MW/ =30.3MW/ 


(b) 


х8.2х ft=0.248 MW/ 


(c) Maximum power 


— 39 МУ/ = 0.364 MW/ 


18. The plate-type fuel element described in Example 6. 1 0 is placed in a test reactor in a region where the 
thermal flux is 5 X 1 013 neutrons/cm2-sec. Calculate: (a) the heat production rate in the fuel; (b) the heat flux 
at the cladding surface; (c) the difference in temperature between the center of the fuel and the cladding surface. 


[sol] 
(a) 
‚ М = 2.86x 


#/ =1.3x Меу/ = 2.06х 


q=  AMeW/sec- 7.8 ABtu /Pr 


(b) 
- =3.25х MeV /secc = 1.65 х Btu / Pr-f 





- =9( —— —=7.8А( —r —_ 0.0004 0 


19. The temperature at the center of the fuel rod described in Problem 8.9, where 4" is the largest, is 1 ,220°F. 
Calculate the temperatures at the fuel-cladding interface at the outer surface of the cladding. 


[sol] 
= 6.663 x Btu / Pr-f , q= = — —= 2.306x 


Btu/Pr 
=1.47 х F-hr/Btu , =3.84x 





і) 9 = => 1220-2.306х х 1.47 х = 881.5 Е 
ii) q= —— 
=> 


20. Derive the general solution to Eq. (8.44). [Hint: Write the first two tenns of the equation as-— ~and then 
perform the integration.] 


[sol] 


2]. Show that Eq. (8.47) can be written as where q' is the heat generated per unit length of a fuel rod. [Note: 
This result shows that the temperature difference across a fuel rod 18 directly proportional to the linear power 
density q' The quantity q' is an important parameter in reactor design.] 


[sol] 


— — m -__—_—— —— — 
— — 


22. The gap between U pellets and the inside of a fuel tube can be taken into account in calculations of heat 
flow out of the fuel by defining a gap conductance by the relation 


where q" is the heat flux and Т is the temperature difference across the gap. Values of of about 1 ,000 


Btu/hr- -F are typical for PWR and BWR fuel. Show that the thermal resistance of a gap is where A is the area 
of the fuel. 


= —where A is the area of the fuel 


[sol] 
а= = = Т1) = — 





23. At one point in the reactor described in Problem 8.8, the heat flux is 280,000 Btulhr-ft2, and the outer 
temperature of the cladding is 563°Е. Compute and plot the temperature distribution in the fuel rod at this 
location (a) ignoring the presence of the gap between the pellets and the cladding; (b) using the gap conductance 
(see Problem 8.22) of 1 ,000 Btu/hr- . 


[sol] 


(a) fuel region 
, 


T(r) = -— 


B.C 
=0 р = q(+) = 2 
= 563 + 2x0.0235x14.583x280000x(6.03x = 4509 
q= 2 


(b) If gap exist 
1) fuel region 
а= ~~ 
2) gap region 





3) cladding region 
T = 0 => Т=1п+ 
B.C i)T( 
= lnr + 
1) ТО = 
-Cln- 





-1139-5.247х x 6.029 х = 822 


24. The fuel rod in the reactor described in Problem 8.8 with the highest heat flux contains uranium enriched to 
2.28 w / o. (a) Calculate the maximum heat-generation rate in this rod when the reactor is operating at design 
power. (b) What is the maximum thermal neutron flux in this rod? [Note: In part (b), use the average fuel 


temperature from Problem 8.23.] 


[sol] 
(а) q= ——— 


b) ( ——— e = 


= 4128.8 -3.19x x x —=-2634.97= 1719 К - 1456 
From problem 8.8 v(u) = 1.65x , (и) = ——— 
N( 


— —=9.954х1 


Btu / Pr-fx 


= 1.15x neutrons/ sec 


25. A slab of iron 5 cm thick, forming part of the thermal shield of an LWR, is exposed to y-rays having an 
average energy of 5 MeV and incident intensity of 1.5 10“ y-rays/cm^-sec. The side on which the radiation 
is incident is held at 585?F; the opposite side is at 505?F. (a) At what rate is energy deposited per ft^ of the slab? 
(b) compute and plot the temperature distribution within the slab. (c) Calculate the heat fluxes at both faces of 
the slab. | Note: For simplicity, assume that the y-ray absorption is exponential with = 0.18cm ,.] 


[sol] 
a) 47-, -E-‘U=eU2*-=,. E-U 
-е 155 х 1.5 x 10'°#/ст?:вес x 3MeV x 0.18cm'! x [(1.548 x 10?Btu/hr-ft)/(MeV/cm^sec)] 


= 5.098 x 107 Btu/hr-ft? 


b) ГТ-0 
d'T/dx?-0 — T-C,X4C, 
B.C  T(0)-To "CT, 


T(x-a)=Ta ~. C;=T(x-a)/a 
^. T=X-( Ta = То)/а+То 
c) 4”--К-К(То-Тх)/а 
From the Table IV.6, . k-10.5/hr-ft-?F 
7. 4”=10.5(585-505)/5-0.03281=5.12 х 103 Btu-ft 


26. Verify that Eq. (8.59) is a solution to Eq. (8.58) and satisfies the given boundary conditions. 


[sol] 

d^T/dx^4—e *-0 

dT/dx- С,+-е *=0 

— T= C,X4C,- -eU* 

B.C i) T(0)-T, Т,-С:- 
ii) Т(а)=Т, Т,-Сүа»С,- -е 

T- D- —Cya+ -е”= (qeUX1).0€,a 


^. Cy=te™-1)( T, - T/a 





С, = T; +— 
^. T=[¢e-1)( T, - T)/a]x + Т, + = a 


27. y-radiation is incident on a slab of thickness a. (a) Derive an expression giving the condition under which the 
temperature in the slab passes through a maximum. (b) In the special case where the two sides of slab are held at 


the same temperature, show that the maximum temperature occurs at 


Provided a 1. 


[sol] 
а) d^ T/dx^^-e "^20 
dT/dx- Cite” _........................ @ 
B.C.dT/dx[x-x,]-0 ^ — C,2--e"*" 
From (2) 
T= С,Х-С,- -eU* 
B.C. 


Т(х-0) = Т, ET C;-— 
“Т--е“"-Х-Т, a-e* 


b) Т,-Т, 

Т(а)= Т,-1, Т, G- uxe ^m, gU) 
(1- ихе". e UX)29 
eU*-1-ua40.5(ua)^-1/6-(ua)? 

1- e*- па-0.5(па)?+1/6-(па)* 

ихе U*™=1- eUX- ца-0.5(иа)?+1/6-(иа)? 
eUX™_]-ya+0.5(ua)-1/6-(ua)” 
-UXm=In[1-ua+0.5(ua)-1/6-(ua)’] 


From In(1-X)-X-X?/242X^/31----- 

In[1-ua40.5(ua)-1/6-(ua)^]- -0.5(ua)+1/6-(ua)” -0.5[-0.5(ua)-1/6:(ua)?? 
-0.5(ua)+1/6-(ua)?-1/8-(ua)” 

-0.5(ua) -1/24-(џа)? 


7. Xm=-1/u[-0.5(ua) -1/24-(иа)?]=а/2(1-иа/12) 


28. The velocity profile 1,000-psia water flowing through an insulated 8-in I.D.(inner diameter) coolant pipe of 
a research reactor can be fit by the expression 


У = 60-5401 
Where T is in °F. (a) Compute the bulk temperature of the water. (b) compute the average velocity of the water 


in the pipe. (c) Compute the value of the Reynolds number. (d) How much water, in Ib/hr, is flowing through the 
pipe? The water density is constant across the pipe at 61.8 lb/ft’ and = 1.270 Ib/hr-ft. 


[sol] 
a) If p, Cp are constans 
Tb= = 
= / (rc-8/2n-ft/121—0.333ft) 
-[155/2r?-245/4(0.333)^--6480(0.333)^][1/2-9/4(0.333)? 
= 128.33°F 
b)V= = 


= =30ft/sec=108000 ft/hr 


c) u=1.270 Ib/hr-ft 
р-61.885/ЛС 


р,=4 =2rc=0.667ft 
Re=—= — — ——23503622 
d) w= p v Ac =r p v re? 
= пх 61.8 x 108000 x (0.333)? 
= 2.33 x 10sIb/hr 





29.compute and plot the temperatures of the water, the outer surface of the cladding, and the fuel center as a 
function of distance along the hottest channel of the reactor described in examples 8.3 and 8.5 under the 
conditions given in Example 8.7. 


[sol] 
H= 12ft 


T f=0.42/2 x 1/12=0.0175ft 


Гс- (0.42+0.024)/2 x (1/2)20.0185ft 
Tm=3970°F 

Tt=543°F 

Kf=1.1 Btu/Hr-ft-°F 

Кс-3.0 Btu/Hr-ft-°F 


Tm-Tf-qR; К,=1/4лНК,) 
Tf-Tc=qR, R,-In(r/rj)/2xHK,) 
Te-Tb=qR,, R,=1/27K,r,H) 


. Tm-Tb=q(R;R.R}) 
К,=6.0286 x 10°°F-hr/Btu 
В.=2.4567 х 10°°F-hr/Btu 

q =4.66 х 107 Btu/hr-ft 

а=лг Н q —7(0.0175)? x 12 x 4.66 x 10725.38 x 10^ Btu/hr 
3970-543-5.38 x 10? x (6.0286 x 10742.4567 x 107+R,,) 
7. R,=9.5618 х 10? ?F-hr/Btu 
hy-1(1/2z x 0.0185 x 12 x 9.5618 x 107)-7.4977 x 103Btu/hr-ft? 


Tc=Tb+qRo=543+5.38 x 10? x 9.5618 х 10° 
=594.4°Е 
Tf=Te+qRc=594.4+5.38 x 10° x 2.4567 x 10°=726.6°F 


30. For the hottest channel in the reactor described in problem 8.9, Calculate the (a) exit temperature of the 
sodium; (b) maximum temperatures of the fuel and the cladding surface. [Note: Take the heat transfer coefficient 
to be 35,000 Btu/hr-ft-^F.] 


[sol] 
Q-mAhzAH-pAcAZCpA Tb 
Q-(dQ/dt)=q= mAh 
= pAc(AZ/At)CpATb= p v Ac Cp ATb 
q-q"AfdZ-q';dzzmCpA Tb “.mCp(dTb/dz)=q’ 


If q'(z)27q'acos(3- , mCp( J= q's4.cos( > 


Ть(®)-Тъһ= (q',,/7mCp)cos | dz 

^Ty(z)2 Ть *(H:q'54,/mCp:-7z)[sinO-1] (m=) 
Ть(®)= Tymt+(H-q’max/ o Cpz)[sinO-1] 

a) Ty Z Ty max=ThO=Tpmt2q” maxVe/ л © Cp 

rj-(0.158/2)-1/12-6.583 x 10° ft 

H=30.5/12=2.542ft 

“Жил ге H= л(6.583 x 10°)* x 2.542=3.46 x 10 ^ft 

г.=(0.158/12)/0.005-(1/12)=7 х 10° 

Ас=(0.21/12)? -x(7 х 107):1.523 x 10^ ft? 


p(Na)=55.06 Ib/ ft 
V=31.2ft/sec 


~ @= p VAc=55.06 x 31.2 x 1.523 x 107=0.2621b/swe=943.21b/hr(per channel) 
Cp=0.3116 Btu/lb-°F from table IV.5 


q” max=6.663 х 10°Btu/hr- t? 

Tpm, =572°Е 

:.Т,-572(2 х 6.66310? x 3.46 x 107) ( x x 943.2 x 0.3116)=1071.4°F 
b) Te max= Tomt Q^ max УгКһЇ(1+ )/]- 

R,=1/hA=1/(35000 x 27 x 7 x 10? x 2.542)=2.556 x 10 ^F-hr/Bt 
а= л © C, Ёу-л x 943.2 x 0.3116 x 2.556 x 107-0.236 


Т. max=572+6.663 х 10° x 3.467 x 2.556 x 107[(14+V1+2.0367)/2.036]=1078.2°F 


Tim тах” Tomt Q^max УФВ (4 

К,=1/4 nH К,=1/(4л х 2.542 x 20)-1.565 х 10°°F-hr/Bt 

R.-In(r. rj)2xH К,= In(7 6.583)2л x 2.542 x 10-3.845 x 10 ^F-hr/Bt 
/.В=В,-+В.+В,=2.205 x 10°°F-hr/Bt, 2x © C, R= m x 943.2 x 0.3116 x 2.205 x 10”-2.036 
"Т, max=572+6.663 х 107 х 3.46 х 107 x 2.205 x 10°[(1+V1+2.0367)/2.036]=1388°F 


31. For the hottest channel in the reactor described in problem 8.7, Calculate, calculate, assuming the reactor is 
bare, the (a) exit temperature of the water; (Б) maximum temperatures of the fuel and the cladding surface. 


| Note: Take h = 7,600 Btu/hr- ft-?F.] 


[sol] 
a) Т, max= Tpn,* 2G) ux У; /ло Ср 


Tom =496°F 9” = 1.3 x 10° Btu/hr-ft 

H=91.9in=7.658ft 

r,=0.298/2 in =1.242 х 10% ,  Vf= n x (1.242 x10’) x 7.658=3.709 x 10°ft° 
o=3.4 x 10°Ib/hr 23142=1469.19 lb/hr (per channel) 

Cp=1.15 Btu/Ib?F 

2) Ty max=496+(2 х 1.3 x 107 x 3.709 x 10°)/( x x 1469.19 x 1.15)=514.2°F 


b) Т, max= Tpm,+ + Ч? пах У: В, [(1+ y] 
R,zI/hA-1/[7600 x 2л х (0.298/2+0.021) x 7.658=1.93 x 10 ^^F-hr/Bt 
a= по C, Ёу-лх 1469.19 x 1.15 x 1.93 x 107-1.025 


Т. max=496+1.3 x 107 x 3.709 x 10? x 1.93 x 1071(1-41-1.0252)1.025|-518.19Е 


с) Tm max— Tpm,+ + Q аа У:К[(1+ y] 
R,=1/4 nH К,=1/4л x 7.658 x 2.5)24.157 x 10°°F-hr/Bt 


К.=1іп(г. гд2лН k;- In{[[(0.298/2)+0.021](0.298/2)]}27 x 7.658 x 10=2.74 x 10°°F-hr/Bt 


K,=2.5 Btu/hr-ft-°F 
K.=10 Btu/hr-ft-°F (for zircaloy) Assuming cladding is zircaloy 


7. R=Rrt+R,.+R),=4.624 x 10°°F-hr/Bt 
=T © C, R= 1 x 1469.19 x 1.15 x 4.624 x 10°=24.54 
. To max=49641.3 x 10’ x 3.709 x 10? x 1.93 x 107[(14+V1424.54")/24.54]=728.2°F 


32. The heat flux at the surface of fuel rods in an experimental reactor varies with position as indicated in Fig. 
9.15. The fuel rods are 0.5 in. in diameter and 6 ft long. Water enters the reactor core at a temperature of 150F 
and flows through the core at a rate of 350 Ib/hr per rod. The heat transfer coefficient 15 1,000 Btu/hr-ft2-F. 
Calculate the temperatures of the coolant, fuel rod surface, and fuel rod center at the entrance of the channel, 3 ft 
up the channel, and at its exit. 








| 4 
2 200,000 | --———m Zz 
Е 
2 
7» 100,000 
= 
0 3 б Figure 8.15 Heat flux v. distance 
Distance along channel, ft for the fuel rod in Prob. 8.32. 
[sol] 


33. Water at about 2,000 psia and 613F leaves a PWR vessel vis a 29-in I.D. primary coolant pipe with an 
average flow velocity of 4935 ft/sec. (a) What is the mass flow rate this pipe in Ib/hr? (b) What is the Reynolds 
number for the water in the pipe? (c) If the water returns to the reactor vessel at the same pressure and velocity 
but at 555F, how large should the return pipe be? 


[sol] 


34. Convert each of the factors in the computation on the Reynolds number given in Example 8.8 under the 
conditions given in that example. 


[sol] 


35. Calculate the mass flow rate in Ib/hr at which the coolant passes up a single channel of the lattice described 
in Example 8.8 under the conditions given in that example. 


[sol] 


36. For the reactor in Problem 8.7, compute the (a) average flow velocity in ft/sec; (b) equivalent channel 
diameter; (c) Reynolds number entering and leaving the hottest channel; (d) heat transfer coefficient at the 
locations in part (c). 


[sol] 


‚ C=0.0311 


37. For the LMFBR core described in Problem 8.9, calculate the; (a) equivalent diameter of a coolant channel; 
(b) total coolant flow area; (c) sodium flow rate in Ib/hr; (d) average sodium exit temperature; (e) Reynolds 
number for the sodium entering and leaving (at the temperature in part [d]) the core; (f) average heat transfer 
coefficient for the sodium entering and leaving the core. 


[sol] 


38. Water enters the core of the BWR described in Problem 8.8 at 525F and flows through the core at a rate of 
106.5* 10? ТЫЛ. The heat production rate in the fuel along the hottest channel can be approximated as 


q" (2) = "n cos(z 2/165) 


where z is in inches. The reactor pressure is 1,020 psia. Compute the (a) maximum-to-average power density in 
the z direction; (b) total flow area in the core; (c) average water velocity near the bottom of the core; (d) 
water/fuel volume ratio; (e) equivalent diameter of a coolant channel; (f) Reynolds number near the entrance to 
the hottest channel; (g) heat transfer coefficient for the convective part of the hottest channel; (h) location of the 
onset of local boiling. 


[sol] 


39. Referring to the PWR described in Examples 8.3 through 8.9: (a) Does boiling of any type occur in this 
reactor? (b) Does bulk boiling occur? (c) Does boiling occur in a channel 30 in from the axis of the reactor? 


[sol] (a) PWR —150bar-2176psia 


o 


— Boiling occurs partially. 
(b) — 
+ (ед. 8.5) 
(еч. 8/9) ° 
(с) — 


40. Starting with the curves of а" апа shown in Fig. 8.12, calculate the curves 


[sol] 


41. (a) Using the results of Problem 8.31 and the Bernath correlation, compute and plot qc" as a function of 
position along the hottest channel of the reactor in Problem 8.7. (b) Compute and plot the DNBR for this 


channel. (c) What is the minimum DNBR? 
[sol] — 


—PWR 150bar 
V=11.97ft/s 


(b) шинж —— 


(c) 


42. Figure 8.16 shows the flow quality as a function of position along the hot channel of the reactor in Fig. 8.12. 
The mass flux for this channel is (a) Compute and plot дс" and ће ОМВК along the 
channel. (b) What is the minimum DNBR? 


[sol] (a) = 


—Fig. 8.16, 


(b) 


43. The maximum heat flux in the HTGR described in Problem 8.6 is 120,000 and the hot channel 
factor is 2.67. (a) What 16 the average heat flux? (b) What is the total heat transfer area? 


[sol] (a) E 


(b) - 


44. The hot channel factor for the LMFBR in Problem 8.9 is 1.85. Calculate the: (a) total heat transfer area in 
the core; (b) average heat flux; (c) maximum heat flux. 


[sol] (a) 


(b) = 


(с) — 


45. The power distribution in the axial direction of the 12ft long fuel rods of a PWR 18 well represented by 


222 


Ч , 


where z is in inches. The minimum DNBR is 2.1 and occurs 90 in up the hottest channel, where the critical heat 
flux is computed to be The hot channel factor for the reactor is 2.78, and the total heat 


999 


transfer area is 48,000 Calculate (a) qma"; (b) дау"; (c) the operating power of the reactor. 


[sol] (a) — 


(b) = 


(c) 


